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I .  SUMMARY 

Theoretical  analyses  were  performed  on  several 
different  types  of  diffusion  flames  to  study  the  flame 
radiation  effect.  In  the  first  problem,  a  soot  formation 
and  oxidation  scheme  was  incorporated  into  a  turbulent 
diffusion  flame  model  adjacent  to  a  solid  fuel.  The 
computed  results  for  the  natural  convective  fire  showed  good 
agreement  with  experimentally  measured  solid  fuel  burning 
rate.  Soot  radiation  increased  its  importance  with  flame 
height.  When  flame  length  is  greater  than  one  meter,  the 
local  radiative  heat  flux  exceeded  that  by  convection.  In 
the  second  problem,  matched  asymptotic  expansions  were 
employed  to  study  the  spherical  diffusion  flame  around  a 
droplet  or  solid  particle  with  flame  radiation.  It  was 
found  that  the  importance  of  radiation  increased  with 
droplet  radius.  The  theory  predicted  that  there  was  a 
maximum  droplet  or  particle  size  above  which  a  spherical 
flame  could  not  be  supported  due  to  radiative  loss.  In  the 
third  problem,  the  thermophoretic  motion  of  small  particles 
(e.g.,  soot)  were  studied  in  a  stagnation-point  laminar 
diffusion  flame.  It  was  found  that  both  the  thermophoretic 
motion  and  the  Brownian  particle  diffusion  can  have  a 
profound  effect  on  the  particle  concentration  distributions 
which  in  turn  will  affect  the  radiative  heat  transfer. 
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II. 


RESEARCH  OBJECTIVE 


The  main  objective  of  this  research  is  to  study 
the  effect  of  flame  radiation  on  diffusion  flames,  in 
particular  those  related  to  solid  and  condensed-phase  fuels. 

III.  RESEARCH  ACCOMPLISHMENTS 

The  major  accomplishments  in  the  research  program 
will  be  briefly  described  in  this  section.  Detailed  results 
are  contained  in  the  Appendices  and  the  list  of  publications 

(A)  The  Effect  of  Soot  Radiation  on  the  Combustion  of  Solid 
Fuels  in  Turbulent  Boundary  Layers 

The  first  part  of  the  effort  dealt  with  the 
adoption  of  low-Reynolds-number  K-e  model  for  solid  fuel 
burning.  In  this  work,  the  chemical  kinetics  are  assumed  to 
be  infinitely  fast  and  the  mixture  fraction  is  assumed  to 
obey  a  beta  probability  density  distribution.  This  work 
shows  that  the  low-Reynolds-number  K-c  approach  is  more 
physically  revealing  than  the  wall  function  method  in 
analyzing  the  surface  boundary  layer  problem.  By  adopting 
the  low-Reynolds-number  model  where  the  same  equations  apply 
all  the  way  to  the  wall,  we  have  discovered  a  number  of 
limitations  and  problems  which  are  hidden  in  the  wall 
function  approach.  To  name  a  few,  the  equation  for  mixture 
fraction  is  no  longer  homogeneous,  as  assumed  in  all  the 
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previous  studies.  The  evaluation  of  mean  mass  fractions  for 
fuel  and  oxidizer  requires  certain  assumption  of  small 
disturbance.  The  boundary  condition  for  g  (square  of 
mixture  fraction  fluctuation)  requires  the  specification  of 
K  (turbulent  and  kinetic  energy)  at  the  wall.  With  a 
pyrolyzing  solid  surface  the  value  of  K  can  be  an  issue. 

W 

Most  details  can  be  found  in  Publication  List  No.  1. 

In  the  second  part  of  the  effort,  finite-rate 
chemistry  was  added  through  an  eddy-dissipation  model.  In 
addition  to  convection  and  diffusive  contribution,  a  soot 
formation  and  oxidation  scheme  is  also  included  in  the  soot 
transport  equation  to  compute  soot  concentration.  Radiative 
heat  transfer  is  calculated  by  using  a  two-flux  method.  The 
model  was  applied  to  the  vertical  burning  of  a  PMMA  slab 
where  there  were  measured  radiation  data. 

In  the  step-by-step  comparison  with  experiment,  it 
is  found  that  low-Reynolds-number  K-c  turbulence  model  can 
produce  good  heat  flux  at  surface  without  any  modification 
in  pure  natural  convection  along  a  heated  plate.  Tie 
modification  of  the  low  Reynolds  number  term  to  accommodate 
surface  mass  addition  is  made  in  this  study.  The  complete 
model,  when  applied  to  the  turbulent  wall  fire,  produces 
good  burning  rate,  radiant  heat  flux  and  average  soot  volume 
fraction.  The  soot  radiation  is  demonstrated  to  be 
important  in  determining  the  burning  rate  of  solid  fuel. 
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More  details  can  be  found  in  Publication  List  No.  2,  which 
is  attached  in  Appendix  A. 

(B)  Structure  and  Extinction  of  Diffusion  Flame  with  Flame 
Radiation  (Publication  List  No.  3,  Appendix  B) 

Using  droplet  combustion  as  a  model  problem,  and 
capitalizing  on  the  temperature-sensitive  nature  of 
radiative  heat  transfer,  the  structure  and  extinction  of 
diffusion  flames  with  flame  radiation  is  studied  via 
multi-scale  activation  energy  asymptotics.  The  flame 
structure  analyzed  consists  of  an  0(c)  reaction  zone 
embedded  within  an  0(5)  radiation  zone  which  in  turn  is 
situated  on  the  0(1)  diffusive-convective  flow  field,  where 
c  «  5  «  1.  The  analysis  yields  the  structure  equation  for 
the  reaction  zone,  which  can  be  cast  in  the  same  form  as 
that  of  Linan's  adiabatic  diffusion  flame  problem  such  that 
his  extinction  results  can  be  readily  used.  Present  results 
show  that  radiative  heat  loss  promotes  flame  extinction  in 
general,  as  expected.  Furthermore,  it  can  also  lead  to  the 
phenomenon  of  dual  extinction  turning  points  in  which  flame 
extinction  due  to  reactant  leakage  and  thereby  kinetic 
limitation  occurs  not  only  for  sufficiently  small  droplets, 
as  is  well  established,  but  also  for  sufficiently  large 
droplets  as  a  result  of  excessive  heat  loss  from  the 
correspondingly  large  flame.  Consequently,  there  exist 
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diffusive-reactive-radiative  systems  for  which  steady 
combustion  is  not  possible  for  all  droplet  sizes.  An 
estimation  of  the  dimensional  radiative  extinction  droplet 
size  is  also  given  for  the  sample  system  studied.  More 
details  can  be  found  in  Publication  List  No.  3  which  is  also 
attached  as  Appendix  B.  This  work  has  been  collaborated 
with  Dr.  B.  H.  Chao  and  Professor  C.  K.  Law  of  Princeton 
University . 

(C)  Concentration  Distribution  of  Particles  in  a 
Thermophoretically  Affected  Flow  Field 

The  magnitude  of  soot  radiative  emission  depends 
not  only  on  the  soot  concentration  level  but  also  depends  on 
where  the  soot  locations  are  in  the  flame  (because  of 
different  temperature) .  Therefore  the  soot  particle  profile 
is  important  for  radiative  calculation.  Since  it  is  known 
that  soot  particles  are  subject  to  thermophoretic  motion  in 
flames,  we  have  undertaken  this  effort  to  investigate  the 
effect  of  thermophoresis  and  particle  Brownian  diffusion  on 
particle  profiles. 

In  the  first  effort,  a  theoretical  analysis  is 
performed  to  study  the  distribution  of  small-size  particles 
in  the  stagnation-point  region  of  a  hot  wall  (no 
combustion; .  Because  of  the  thermophoretic  motion,  a 
particle  stagnation  point  is  created  in  the  interior  of  the 
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flow  which  can  produce  singularity  when  particle  diffusion 
is  neglected.  By  including  Brownian  diffusion  and  using 
matched  asymptotic  expansions,  the  structure  of  the  inner 
diffusion  layer  near  the  particle  stagnation  point  is 
resolved.  The  results  presented  in  Appendix  C  illustrate 
how  the  diffusion  layer  is  evolved  and  its  structural 
variation  with  a  number  of  parameters. 

In  the  second  part  of  the  effort,  a 
two-dimensional  counterflow  diffusion  flame  stabilized  in 
the  forward  stagnation  region  of  porous  cylinder  is 
analyzed.  Fuel  is  ejected  from  the  cylinder  with  a 
specified  blowing  velocity  and  diffuses  into  an  oxidizer 
flow  of  specified  stagnation  velocity  gradient.  In  this 
configuration,  particles  (e.g.,  soot)  of  an  assumed  average 
dimension  form  at  a  global  rate  depending  on  the  local  fuel 
concentration  and  temperature.  Thermophoresis  gives  rise  to 
many  different  flow  patterns  for  the  particles  depending  on 
the  flame  structure  and  the  magnitude  of  thermophoretic 
coefficient.  In  some  cases,  multiple  stagnation  points  for 
particle  flow  exist. 

Brownian  diffusion  needs  to  be  included  in  this 
theoretical  treatment  in  order  to  avoid  singular  behavior 
(infinite  particie  concentration)  at  the  stagnation  point 
under  certain  conditions.  Because  the  Schmidt  number 
associated  with  the  Brownian  diffusion  is  very  large  a 
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special  numerical  scheme  was  developed  to  solve  the  mass 
conservation  equation  for  the  particles.  The  numerical 
scheme  was  implemented  in  two  stages  and  allowed  to 
integrate  the  equation  even  in  situations  for  multiple 
turning  points  (points  where  particle  velocity  changes 
direction) .  This  was  achieved  by  making  use  of  an  adaptive 
mesh  refinement  scheme. 

For  cases  along  the  center  streamline  where  the 
particles  are  convected  toward  the  stagnation  point,  a 
boundary  layer  behavior  is  observed.  When  the  product  of 
Prandtl  number  a  the  thermophoret ic  coefficient  is  greater 
than  a  critical  value  (nearly  unity  -  slightly  influenced  by 
chemical  kinetics) ,  a  local  peak  of  particle  concentration 
is  predicted.  For  cases  where  particle  velocity  is  from  the 
fuel-rich  side  toward  the  flame,  the  particles,  generated  in 
the  fuel-rich  side  (such  as  soot)  ,  can  leak  through  the 
flame  and  appear  on  the  fuel-lean  side. This  can  occur  only 
when  the  stoichiometric  fuel-oxidizer  mass  ratio  is  less 
than  unity.  The  detailed  profile  of  particle  concentration 
in  counterflow  diffusion  flames,  therefore,  depends  on  many 
parameters.  More  details  can  be  found  in  Publication  List 
No.  5. 
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ABSTRACT 

A  comprehensive  model  for  turbulent  boundary  layer  combustion  adjacent  to  solid  fuels  is  estab¬ 
lished  and  employed  to  study  the  influence  of  thermal  radiation  on  the  burning  rate  of  solid  fuels.  A  low- 
Reynolds-number  two-equation  k  -  e  model  is  used  for  flow  transport.  The  treatment  of  turbulence  and 
chemistry  interaction  follows  the  eddy-dissipation  model.  Radiative  heat  transfer  is  calculated  by  using  a 
two-flux  method,  (n  addition  to  convection  and  diffusive  contribution,  the  soot  formation  and  oxidation 
are  also  included  in  the  soot  transport  equation  to  compute  the  soot  concentration. 

Comparisons  between  model  calculations  and  experimental  data  are  carefully  carried  out  for  a  verti¬ 
cally  burning  PMMA  fire.  Consistent  with  the  experiment  results,  the  model  calculation  shows  that  radia¬ 
tion  to  solid  fuel  is  shown  to  be  the  most  important  factor  in  determining  the  solid  pyrolysis  rate. 
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INTRODUCTION 

In  t lie  past  decade,  a  large  amount  of  attention 
was  given  to  fire  study  Detailed  experiments  1 22.  2t>.bj 
detnunstrated  that  although  convection  is  the  dominant 
iteat  transfer  mode  in  small  scale  fires,  radiation  heat 
iranster  is  the  dominant  mode  in  large  fires.  For  the 
burning  ot  hydrocarbons  or  polymers.  it  has  long  been 
recognized  that  soot  particles  are  the  important 
radiation  sources  in  the  flame.  Recently,  it  was  pointed 
out  (Id.  24.  5)  that  fuel  vapors  or  their  decomposition 
products  also  play  important  roles  in  the  flame  radia¬ 
tion.  Therefore,  in  order  to  model  the  burning  behavior 
of  large  combustion  systems  properly,  a  theoretical 
model  should  be  able  to  determine  the  amount  of  soot, 
tuei  vapor  and  combustion  product,  their  distribution 
in  the  flume,  fame  shape  and  sue. 

Theie  are  many  theoretical  efforts  to  model  the 
solid  fuel  burning  rate  when  radiation  is  significant.  Al¬ 
though  these  radiation  models  vary  in  their  complexities, 
they  are  all  based  on  very  drastic  assumptions.  For 
example,  in  studying  turbulent  wall  tires,  two  radiation 
models  are  assumed  by  Tamamni  [ 3 2 1  .  In  the  first 
model,  radiation  power  is  assumed  to  be  a  constant 
traction  of  the  energy  release  rate.  In  the  second  model, 
the  radiation  is  postulated  to  emit  from  a  thin,  constant- 
temperature  layer  of  particles  at  the  flame  front  and  the 
value  of  ennssivity  of  the  flame  sheet  is  chosen  empiri¬ 
cally.  In  the  work  of  Methochianakis  [23 1 .  the  flame 
ernissivity  is  assumed  to  be  related  to  the  unburned  fuel 
in  an  exponential  manner.  In  Ref.  11,  radiative  heat 
tlux  is  simply  taken  as  constant  mean  value. 

The  purpose  of  this  article  is  to  present  a  more 
comprehensive  model,  compare  its  calculation  with 
available  data  and  evaluate  the  importance  of  soot 
radiation.  This  model  includes  several  physical  sub¬ 
models  which  describe  the  turbulence  phenomena, 
turbulence-chemistry  interaction,  radiation  heat  transfer 
and  the  soot  formation  and  oxidation. 

MATHEMATICAL  MODEL 

1 .  Mean  Transport  Equation 

The  physical  situation  studied  is  schematically 
shown  in  Fig.  1.  A  turbulent  buoyant  flame  is  sustained 
by  the  heat  feedback  to  the  solid  fuel  due  to  convective 
and  radiative  heat  fluxes  from  the  high  temperature 
gases  in  the  (lame.  Some  assumptions  are  made  for  easy 
treatment  of  the  problem:  (a)  the  average  flow  is  two- 
dimensional  and  steady,  (b)  the  binary  diffusion  coef¬ 
ficients  for  all  pairs  of  species  except  soot  and  soot 
precursor  are  equal,  (c)  mass  diffusion  due  to  pressure 
gradients,  temperature  gradient  are  neglected,  (d)  the 
soot  and  soot  precursor  are  treated  as  minor  species, 
their  heat  of  formation,  combustion  heat  and  oxygen 


Fig.  1.  Schematic  of  the  problem  studied 

consumption  are  neglected,  (e)  the  thermophoretic 
motion  of  soot  is  neglected,  (0  the  specific  heat  Cp 
is  constant,  (g)  Lewis  numbers  except  for  soot  and  soot 
precursor  are  equal  to  one.  (h)  the  ideal  gas  law  applies 
to  the  gas  mixture. 

Following  the  mass-averaging  procedure  and  order 
of  magnitude  analysis,  the  following  two-dimensional 
boundary-layer  equations  are  obtained 


b(pu) 

bx 


biH  v) 
by 


=  0. 


(It 


bipuu)  b{J5vu)  b  ,  bu  , 

ax  ~ y  by  'M'//  av  '  +  ^  ~  Paa)' 

,3, 

bx  by  by  Sc,  Scj,  by  '• 

b(puh, )  b(pvhr)  _  b  .  p  u,  bh,  . 

bx  by  by  Pr  Prt  by 

Pr  Pr / 


P  = 


by 

t  pixu  +  SR  . 
pRT 


by 


nr 


(4) 

(5) 


where  k  is  the  turbulent  kinetic  energy  and  SR  is  the 
heat  addition  due  to  radiation  heat  transfer.  All  the 
other  symbols  are  described  in  the  Nomenclature. 
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2.  Turbulence  Model 

Two-equation  k  -  e  model  is  u^ed  to  calculate  the 
turbulent  viscosity  and  thus  close  the  averaged  equa¬ 
tions  Instead  of  using  wall-functions  { 1 5  j  to  estimate 
the  transfer  ot  momentum,  mass  and  energy  at  solid 
surface.  Clnen  s  low-Reynolds-nutnber  k  -  e  model  [Sj 
is  employed  to  make  the  calculation  right  up  to  the  solid 
nurtace  The  transport  equations  tor  turbulent  kinetic 
energy  \ .  and  the  dissipation  rate  e  are. 
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P,  =  0.09  [I  -  t'.xp  (  -  C3  l'+  ) 

/,  -  I  35  f 
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The  real  turbulent  kinetic  energy  dissipation  rate 
e  is  related  to  i  by  e  =  i  +  2 vkjy1 .  In  this  study,  the 
constant  C3  has  been  modified  from  its  original  value 
0.0 1 15  (5 1  to  account  for  the  effect  of  mass  transpira¬ 
tion  through  the  solid  surface  as 

C3  =  0.01 15/(1  +  Qcvwju*)  (8) 

This  point  will  be  discussed  in  the  later  section. 

The  buoyancy  generation  terms  in  the  equations 
(6)  and  (7)  are  neglected  Lockwood  et  al.  [18j 
compared  two  conventional  k  -  e  turbulence  models 
with  and  without  buoyancy  generation  terms,  they 
found  that,  for  turbulent  round  diffusion  (lame,  the 
results  are  insensitive  to  these  terms.  A  application  ol 
equations  (6)  and  (7)  to  natural  convection  from  a 
vertical  isothermal  surface  also  supports  this  finding. 
The  details  are  presented  in  the  Results  and  Discussion 
Section. 

3.  Combustion  Model 

The  combustion  of  fuel  is  treated  as  a  single  step 


irreversible  reaction: 

I  Kg  [tue!|  +  s  kg  [ 0 J  -»  [  I  +  s  |  kg  [product)  [Oi 

The  nonlinear  interaction  between  the  turbulence 
and  reaction  makes  the  time  averaging  of  the  reaction 
rate  very  difficult  to  evaluate.  For  high  local  Reynolds- 
number  U‘,Ve)  region,  we  assume  that  the  local  rate  ot 
reaction  is  limited  by  the  molecular  mixing  between 
reactants  and  tollow  the  approach  of  Magnussen  |21] 
by  expressing  the  reaction  rate  as: 

W,  =  :3bP  4  —min.  |  V,  .  — —I- T  1—  <  I0t 
k  k  !  s  I  •  7* 

where 


For  low  local  Rey  nolds-number  region,  we  assume 
the  reaction  rate  is  chemically  controlled  and  is  ex¬ 
pressed  as  (9) : 

i mi 

The  criterion  of  low  Reynolds-number  is  quite 
arbitrary,  Eq.  ( 1 1 )  is  used  when  local  Reynolds  number 
is  less  than  50. 

4.  Soot  Model 

The  processes  governing  the  rates  of  formation  and 
oxidation  of  soot  are  quite  complex  and  quantitative 
models  to  describe  these  processes  are  not  yet  fully 
developed.  Quasi-global  models  [34,  14.  8|  have  been 
developed  to  describe  the  process  of  soot  formation  via 
a  few  overall  steps.  Tesner  et  al.  [34]  proposed  that 
soot  formation  occurs  in  two  steps.  The  first  stage  re¬ 
presents  formation  of  radical  nuclei  (soot  precursors), 
and  the  second  stage  the  soot  particle  formation  from 
these  nuclei.  A  modified  Tesner's  model  [10)  is 
adopted  in  this  work.  The  model  describes  the  forma¬ 
tion  rate  of  soot  precursor  as 

~TT~  A  Yftxp  (~£niRT)  +  (f-  go  «(.V  +  n ) 

1  fo 

(12) 


and  the  formation  rate  of  soot  as 


dN  - 

~  g0  n{N  +  n)  j  K 


(13) 


The  oxidation  of  soot  is  better  understood  and  has 
been  modelled  by  several  investigators  (16,  25.  19). 
It  was  suggested  that  the  soot  oxidation  is  chemically 
controlled  in  many  situations.  However  Magnussen  (20) 
and  Abbas  et  al.  (1 1  argue  that  the  rate  of  burning-out 
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may  bo  controlled  by  turbulent  mixing  exspecially  in 
high  temperature  region.  In  this  study,  we  follow  Abba:, 
el  al.  and  express  the  burning  rate  as  the  smaller  of 
kmctically  controlled  burning  rate  Wi  k  and  mixing 

controlled  burning  rate  IV.  trj. 

tv^  ^  -  min.  | IV.  ...  IV.  ni  ]  (14) 


where 


IV 


/J", 

.1  ~ •  j  x  exp  (  RT\. 


(15) 


tv,  r, 

tv  , ,  = - - -  (Its) 

1  '■  y. 

The  oxidation  ot  precursors  is  assumed,  following 
Magnussens  work,  as 

-r  tv. or 


where  ft'  is  the  nucleus  particles  per  unit  mass  of  mix¬ 
ture. 

5  Radiation  Heat  Transfer 

The  radiation  heat  transfer  is  calculated  by  using 
two-flux  radiation  model  ( 30 j .  The  equation  for  the 
radiative  flux  away  from  the  solid  surface  q  and  the 
radiative  tlux  to  the  solid  surface  q~  are  derived  by 
using  Schuster-Schwarzschid  method  [31  J.  The 
resulting  equations  are: 


bq+ 

dy 


-  a'q +  +  a'F. 


(18) 


a'q  +  a'F, 


(19) 


where  a'  is  the  absorption  constant  and  E-  oT* .  De¬ 
fining  F  -  0.5(rjr+  +  q~)  and  q  =  q+  -  q~ ,  from  Eq. 
(18)  and  Eq.  (19)  we  have 


_d_ 

by 


(~ 
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bF_ 

bv 


)  + 


a'  (E  -  F)  -  0. 


(20) 


2  bF 
^  a’  by 


(21) 


The  radiative  energy  contribution  to  the  source 
term  in  the  total  enthalpy  equation  is  expressed  as 

SR  -*(F-E)  (22) 

The  absorption  constant  a'  =  2 Ka  and  the  absorp¬ 
tion  coefficient  of  the  mixture  Ka  is 


Kj  =  I  MJ8  f/v  +  20  pY,  +(6  2  0  004367'  + 

0  898  x  I  O'6  7:  )V_  1 23 ) 

r 

where  J\.  is  file  volume  fraction  of  soot.  The  contribu¬ 
tion  of  -.out  to  Ka  is  taken  from  Ret.  I7.  I  he  contribu¬ 
tion  of"  product  is  deduced  from  data  of  Rapanotu  |29| . 
and  contribution  of  fuel  vapor  is  consistent  with 
available  data  |3) 

6.  Boundary  Conditions 

In  addition  to  the  equations  reported  in  the 
pievious  sections,  the  specification  ot  boundary  condi¬ 
tions  is  required  to  complete  the  mathematical  formula- 
tion. 

At  the  edge  of  the  boundary  layer,  we  have. 


u  =  0 

t24tf) 

y,  =  o 

(246) 
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The  boundary  conditions  at  solid  surface  are: 


3  =  0. 
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The  condition  (  2 5 1 ")  is  exact  only  lor  rigid  wall,  tor 
sjpori/ing  wall,  this  implies  negligible  burning  rate 
fluctuation  Conditions  (25b)  -  (25d)  are  obtained 
from  the  species  mass  balance  at  the  solid  fuel  surface. 
In  condition  (25i).  it  is  assumed  that  the  fuel  regression 
rate  can  be  represented  in  an  Arrhenius  form.  Condi¬ 
tion  (25hl  tollows  from  the  relation  q+  =  awoT^.  +  ( 1  - 
au  )q  '  at  wall  and  finally  the  condition  (25j)  is  obtained 
from  the  energy  balance  at  the  solid  surface. 

Numerical  Method 

The  Patankar-Spalding  finite  difference  procedure 
(2~|is  chosen  for  the  numerical  solution  of  the  system 
ol  equations.  One  hundred  and  fifty  cross  stream  grid 
nodes  are  used  to  obtain  grid-independent  solution 
Over  sixty  points  are  located  within  y+  <  20.  The  tlow 
is  assumed  to  start  out  laminar  and  the  starting  profiles 
are  the  similarity  solutions  corresponding  to  the  pres¬ 
cribed  burning  rate  and  the  radiation  tlux  for  constant 
wall  temperature  situation.  The  turbulent  calculation 
is  initiated  by  introducing  k  and  e  at  a  distance  0.04  m 
from  the  leading  edge.  The  initiating  profile  of  k  and  e 
chosen  are  (22 1 


X  =  Cjr  . 

(26) 

e  =  Cek^-  l 

(27) 

where  Ck  and  Cf  are  constant  and 

/  =  y  for  y  <  0.0008m, 

/  =  0.0008m  for  y  >  0.0008m 

RESULTS  AND  DISCUSSIONS 

The  numerical  calculation  is  conducted  for  PMMA 
fuel.  The  physical  properties  and  model  constants  used 
in  this  work  are  listed  in  TABLE  1.  The  viscosity  p  = 
1 .4573  x  10'7  x  1 10.4  +  T)  varies  with  tempera¬ 
ture.  The  constant  a  in  Eq.  (8)  for  the  low-Reynolds- 
number  k  -  e  turbulence  model  is  set  to  1.75  in  order 
to  get  better  performance  if  there  is  transpiration  at  the 
surface.  A  comparison  was  made  between  model 
calculations  and  the  experimental  data  (13  j  in  the 
forced  incompressible  boundary-layer  flow.  The  results 
are  shown  in  Fig.  2.  In  the  figure,  G  =  is 

defined  as  the  mass  flux  transpired  at  the  solid  surface 
to  the  mass  tlux  at  the  free  stream.  The  model  per¬ 
formance  is  good  near  the  wall.  Since  the  velocity- 
profile  at  large  y+  is  more  sensitive  to  external  parame¬ 
ters,  particularly  the  pressure  gradient.  The  deviation 
at  large  y +  for  G  =  0.0038  may  be  due  to  the  different 
pressure  gradients  used  in  the  numerical  caluclation  and 


Table  I 

Values  of  Properties  and  Model  Constants 


Properties 
or  Model 

Constants 

Value 

Lnits 

CP 

1254 

J  Kg-7k 

Cs 

1463 

J  kg-'k 

Tso 

200 

°k 

Too 

200 

’k 

L, 

1 .05  x  106 

J  kg-'k 

0.48 

s 

1.92 

4 

2.5957  x  10' 

J  kg-:k 

V  0 ,  00 

0  2215 

Ds 

250 

A° 

^soot 

2000 

kg,  m3 

Poo 

1 

atm 

A 

1.2  x  I041 

part/ m3  s 

l'-g 

100 

s'1 

Y,b 

0.1 

K 

2.5 

go 

1  x  10'15 

m3  /s/  part 

En 

180 

kcal/mole-°K 

Sc, 

10 

Scn 

10 

Scj  (i  *  s,m) 

0.7 

Pr 

0.7 

Pr, 

0.9 

Sc, 

0.9 

Ag 

5.047  x  107 

s’1 

Eg 

27 

kcal/mole-°K 

Ew 

30 

kcal/mole-°K 

B 

5.144  x  10® 

Kg/m2  -see 

Ck 

0.02 

8x 

-9.8 

m/s2 

a 

1.75 

A* 

1.302  x  I0‘° 

Kg-°K0S/m3- 

E, 

39.3 

kcal/mole-°K 

in  the  experimental  work.  The  low-Reynolds-number 
turbulence  model  was  also  tested  in  pure  natural  convec¬ 
tion  situation.  In  the  case  of  buoyant  How  near  a 
vertical  heated  plate,  the  resulting  surface  heat  tlux  and 
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Fig.  3.  Local  Nusselt  number  versus  Grashof  number 

the  temperature  profile  are  shown  in  Fig.  3  and  Fig.  4 
respectively.  The  results  are  in  good  agreement  with  the 
experimental  data  and  the  calculation  by  Plumb  and 
Kennedy  (28)  using  k  -  e  -  T'2  turbulence  model 
which  is  also  a  low-Reynolds-number  turbulence  model 
but  with  buoyant  generation  terms  in  x  and  e  equations. 
This  good  agreement  supports  the  finding  of  Lockwood 
[  1 8 1  that  the  buoyant  generation  terms  can  be  neglected 
for  boundary  layer  flows  with  gravity  vector  nearly 
perpendicular  to  temperature  gradient.  Fig.  3  also 
indicated  that  disturbance  introduced  at  different 
Grashoff  number  (2  x  107  and  4  x  108)  leads  to  the 
same  result  in  fully  turbulent  region. 

The  complete  model  was  applied  to  the  turbulent 
buoyant  diffusion  flame  on  a  vertical  surface.  Using  the 
soot  model  constants  (A  =  6.25  x  1040,  Yja  =  0.2  and 
K  =  5),  suggested  by  the  Galant  et  al.  [10) ,  it  is  found  in 
Fig.  5  that  the  discrepancy  between  the  predicted  and 
experimental  burning  rate  [26]  is  large  and  the  cal¬ 
culated  averaged  volume  fraction  (~0.85  x  10'7)  is 
much  lower  than  the  reported  data  (~  0.3  x  1  O'6 )  [6|. 
This  discrepancy  may  due  to  the  fact  that  different  fuel 
is  burned.  The  constants  suggested  by  Galant  et  a),  is 
for  methane  flame,  while  the  experimental  data  is 
obtained  by  burning  PMMA.  Using  the  constants  shown 


Height  (m) 


Fig.  5.  Pyrolysis  rate  as  a  function  of  height 

in  Table  I.  we  could  make  the  calculated  burning  rate 
quite  close  to  the  experimental  data  and  the  calculated 
average  soot  volume  fraction  (  ~  0.27  x  1 0’6)  agrees 
with  the  reported  data. 

Fig.  6  shows  the  ratio  of  radiant  heat  flux  to  the 
surface  to  the  total  heat  flux.  It  indicates  that  the  radia¬ 
tion  heat  transfer  is  a  dominant  heat  transfer  mode  in  a 
big  fire.  The  calculated  radiant  heat  flux  and  convective 
heat  flux  to  the  surface  are  3iso  compared  with  the 
results  of  Orloff  et  al.  [26)  in  Fig.  7.  Orloff  et  al. 
deduced  the  radiant  heat  flux  and  the  convective  heat 
flux  from  experimental  data  using  optically  thin  flame 
analysis.  The  radiant  heat  flux  using  the  gray  flame 
analysis  is  3lso  presented  in  Fig.  7.  From  Fig.  7,  we  can 
see  that  the  agreement  between  the  model  calculation 
and  the  results  of  Orloff  et  al.  is  good.  The  calculated 
convective  heat  flux  decreases  slightly  with  the  height, 
while  the  radiant  heat  flux  increases  with  the  height. 
The  decrease  of  the  convective  heat  flux  is  due  to  the 
decrease  of  gas  temperature.  As  will  be  shown  in 
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Fig.  7.  Radiant  and  convective  heat  transfer  to  the 
surface  versus  height 

parametric  studies,  the  convective  heat  flux  increases 
slightly  m  no-gas-phase-radiation  case  even  the  flame 
temperature  decreases  slightly  with  the  height.  Radia¬ 
tion  in  gas  phase  makes  the  gas  temperature  decreases 
more  and  it  is  this  further  decrease  of  gas  temperature 
that  makes  the  convective  heat  flux  decrease  instead 
tncrease  slightly  with  the  height.  The  increase  of  radiant 
heat  is  mainly  due  to  the  fact  that  the  soot  volume 
fraction  increases  with  height. 

Fig.  8  and  Fig.  9  show  the  typical  profiles  of 
velocity,  temperature,  fuel  mass  fraction,  oxygen  mass 
fraction,  soot  mass  fraction,  turbulent  kinetic  energy 
and  dissipation  rate.  The  soot  mass  fraction  is  much 
larger  on  the  fuel-rich  side  than  on  the  fuel-lean  side.  It 
is  consistent  with  the  experimental  finding  that  the  soot 
layer  is  on  the  fuel  side  of  the  flame  ^one.  Because  of 
the  shape  of  velocity  profile,  the  turbulent  kinetic 
energy  and  the  turbulent  energy  dissipation  profiles 
have  dents  near  the  peak  of  the  velocity  profile.  As 
shown  in  Fig.  10.  further  downstream,  the  turbulent 
kinetic  energy  profile  shows  two  local  maximum  values. 


Fig.  S  Profiles  of  velocity,  temperature,  fuel  mass 
traction  and  oxygen  mass  iiuvtion 
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Fig.  9  Profiles  of  soot,  turbulent  kinetic  energy  and 
turbulent  kinetic  energy  dissipation  rate 
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Fig.  10.  Profiles  of  soot,  turbulent  kinetic  energy  and 
turbulent  kinetic  energy  dissipation  rate 

Parametric  studies  are  also  presented  in  Fig.  11. 
Two  additional  cases  are  considered:  one  without  fuel 
contribution  in  gas-phase  radiation,  and  the  other 
without  gas-phase  radiation  at  all  but  still  with  surface 
radiation.  We  note  that  the  radiant  heat  transfer  from 
soot  is  large  compared  with  that  from  fuel  vapor.  The 
figure  also  indicates  that  in  the  no-gas-phase-radiation 
case,  the  burning  rate  increases  slightly  with  height.  In 
this  case  the  heat  flux  to  the  surface  is  totally  due  to  the 
convection,  so  the  convective  heat  flux  increases  with 
height.  The  reason  why  the  convective  heat  flux  in¬ 
creases  with  height  is  that  the  gas  velocity,  the  boundary 
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layer  thickness  and  therefore  the  tuihuient  Jin'usiv u\ 
mcieasc  with  heigh'.  It  is  not  due  to  the  tlaine  tempera¬ 
ture  but  tite  temperature  profile,  because  the  :na\imum 
gas  temperature  decreases  slightly  with  height. 

In  numerical  calculation,  it  is  found  that  the 
transition  from  laminar  region  to  fully  turbulent  region 
is  quite  sensitive  to  the  triggering  method.  Fig.  Id  shows 
some  results.  In  all  calculations  the  initial  tuihuient 
kinetic  energy  is  kept  the  same  but  the  initial  dissipation 
rate  varies.  The  ligure  indicates  that  depending  on  the 
initial  disturbance,  the  laminar-turbulent  transition  point 
varies.  We  note  that  the  curves  are  nearly  parallel  to 
each  other  at  fully  turbulent  region  and  shift  from  each 
other.  This  shift  is  due  to  the  variation  of  the  laminar- 
turbulent  transition  point  Partial  support  for  this  is 
that  in  fully  developed  region,  at  same  local  burning 
rate,  they  (except  C(  -  0.3  which  is  still  in  early 
turbulent  region)  all  nearly  have  the  same  boundary 
layer  thickness.  The  possibility  of  using  low-Reynoids- 
number  k  -  e  turbulence  model  to  calculate  the  transi¬ 
tion  region  and  also  the  influence  of  the  initial  dis¬ 
turbance  is  still  under  study.  At  the  present  stage,  the 
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Fig.  12.  Pyrolysis  rate  as  a  function  of  height 


transition  portion  only  provides  a  means  >1  leading  to 
Mi  luily  developed  turbulent  region  winch  is  ol  interest. 

CONCLUSIONS 

I  he  pyrolysis  and  combustion  ol  solid  tuel  a 
tuihuient  boundary  layer  jre  determined  by  the 
complex  interaction  ol  various  transport  and  kinetic 
processes  In  otdei  to  assess  the  relative  importance  of 
each  component  process,  a  comprehensive  model  winch 
consists  ol  several  sub-models  is  established  m  this  work 
Although  simple  expressions  are  given  to  the  sub-model 
processes,  the  combined  model  is  juiie  eoiimlex 
Because  the  only  existing  experiment  on  >oiid  'ue! 
burning  in  turbulent  flame  with  detailed  measurement 
of  radiative  properties  is  in  njiural  convection.  :iie 
model  is  ;iisi  validaied  m  this  situation. 

In  the  step-by-step  comparison  with  expenmetu. 
ii  is  tumid  that  C'hien's  iow-Rey nolds-nimiber  x  e 
tuibulence  mode!  can  produce  good  heat  llux  at  surface 
without  any  modification  ;n  pure  natuial  convection 
jlong  healed  plate.  The  moditication  of  the  low 
Reynolds  numbei  term  to  accommodate  surface  mass 
addition  is  made  m  tins  study.  Tiie  complete  model, 
when  applied  to  the  tmbulent  wall  tire,  produces  good 
burning  rale,  radiant  heat  tlux  and  average  soot  volume 
traction.  The  soot  radiation  is  demonstrated  to  be 
important  in  determining  the  burning  rate  of  solid  fuel. 

The  prediction  of  dit fusion  flame  in  a  turbulent 
boundary  layer  depends  on  the  upstream  conditions 
which  are  usually  difficult  to  specify  completely.  In  the 
present  study,  a  laminar-to-turbulent  transition  pro¬ 
cedure  is  adopted.  Bui  [he  resuits  also  depend  on  how 
the  triggering  disturbance  is  specified.  The  determina¬ 
tion  of  the  most  proper  triggering  disturbance  in  t he 
low-Reynolds-number  k  -  e  turbulence  model  to  bridge 
the  laminar  and  turbulent  regions  will  need  more  work. 
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NOMENCLATURE 

A  Model  constant 

Ag  Preexponential  constant 

As  Preexponential  constant 

a'  Absorption  constant 

au.  Lmissivity  of  solid  luel 

B  Preexponential  constant 

Ck  Coetficient 
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Specific  heat  of  the  mixture 

Too 

Ambient  temperature 

Specific  heat  of  the  solid  fuel 

Pso 

Bulk  temperature  of  the  solid  fuel 

c 

«■ 

Coefficient 

t 

Time 

D 

Species  Jiffusivitv 

it 

Velocity  in  the  v  direction 

Soot  particle  diameter 

u* 

E 

=  aT4 

it* 

=  u  it  * 

Activation  energy 

V 

Velocity  in  they  direction 

Activation  energy 

W 

Mean  molecular  weight 

Activation  energy 

Vv, 

Neat  rate  ot  formation  of  species  /  per  unit 

Activation  energy 

volume 

E 

=  0.5  U/  +  +  q~ ) 

*n.c 

Rate  ot  consumption  of  nuclei  per  unit 

J 

Model  constant 

Kc 

Rate  of  consumption  of  soot  per  unit  volume 

f\ 

Volume  traction  of  soot 

Wf.k 

Rate  of  consumption  of  soot  per  unit  volume 

G 

=  wu  pxi'oo 

— 

controlled  by  chemical  kinetics 

(Jr 

Local  Grashof  number 

Km 

Rate  of  comsumption  of  soot  per  unit  volume 
controlled  by  mixing 

g 

Model  constant 

V 

Coordinate  along  the  wall 

g  0 

Model  constant 

Y,o 

Model  constant 

?x 

Acceleration  of  gravity  in  x  direction 

y, 

Mass  fraction  of  species  i 

\ 

Total  enthalpy 

V 

coordinate  normal  to  the  wall 

K 

Model  constant 

+ 

V 

=  yu*jv 

Ka 

Absorption  coefficient  of  the  mixture 

a 

Model  constant 

k 

Turbulent  kinetic  energy 

T* 

Mass  fraction  occupied  by  tine  structures 

Ls 

Latent  heat  of  the  solid  fuel 

e 

Rate  of  dissipation  of  turbulent  kinetic  energy 

l 

Length 

e 

=  e  -  2kv/v‘ 

m 

Mass  flux 

\ 

Thermal  conductivity 

y 

Concentration  of  soot  particle  (part/m3) 

M 

Molecular  dynamic  viscosity 

n 

Nucleus  concentration  (part/m3 ) 

P, 

Turbulent  viscosity 

n' 

=  nip 

P'ff 

=  H  +  P, 

P 

pressure 

V 

Molecular  kinetic  viscosity 

Pr 

Prandtl  number 

p 

Density 

PG 

Turbulent  Prandtl  number 

0 

Planck  constant 

q 

Heat  of  combustion  per  unit  of  mass  of  fuel 

Wall  shear  stres 

q 

<7+ 

=  q+  -  q~ 

Radiative  heat  flux  in  positive  y  direction 

Superscripts 

q~ 

Radiative  heat  flux  in  negative  y  direction 

- 

Conventional  time  mean  value 

R 

Universal  gas  constant 

- 

Favre  mean  value 

s 

5c, 

Stoichiometric  oxidizer/fuel  mass  ratio 

Schmidt  number  of  species  i 

Subscripts 

Sci  t 

Turbulent  Schmidt  number  of  species  i 

f 

fuel 

Sr 

Heat  addition  due  to  radiation  heat  transfer 

i 

species  i 

T 

Temperature 

07 

Oxygen 

Tw 

Surface  temperature 

n 

Nucleus 

A10 
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V  product 

1  vOOt 

r  turbulence 

«•  Wall 

»  Vmbience 
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ABSTRACT 

Using  droplet  combustion  as  a  model  problem,  and  capitalizing  on  the  temperature-sensitive 
nature  of  radiative  heat  transfer,  the  structure  and  extinction  of  diffusion  flames  with  flame  radiation 
is  studied  via  multi-scale  activation  energy  asymptotics.  The  flame  structure  analyzed  consists  of  an 
O(e)  reaction  zone  embedded  within  an  0(5)  radiation  zone  which  in  turn  is  situated  in  the  0(D 
diffusive-convective  flow  field,  where  e«5«l.  The  analysis  yields  the  structure  equation  for  the 
reaction  zone,  which  can  be  cast  in  the  same  form  as  that  of  Lilian's  adiabatic  diffusion  flame  problem 
such  that  his  extinction  results  can  be  readily  used.  Present  results  show  that  radiative  heat  loss 
promotes  flame  extinction  in  general,  as  expected.  Furthermore,  it  can  also  lead  to  the  phenomenon  of 
dual  extinction  turning  points  in  which  flame  extinction  due  to  reactant  leakage  and  thereby  kinetic 
limitation  occurs  not  only  for  sufficiently  small  droplets,  as  is  well  established,  but  also  for  sufficiently 
large  droplets  as  a  result  of  excessive  heat  loss  from  the  correspondingly  large  flame.  Consequently 
there  exist  diffusive-reactive-radiative  systems  for  which  steady  combustion  is  not  possible  for  all 
droplet  sizes.  An  estimation  of  the  dimensional  radiative  extinction  droplet  size  is  also  given  for  the 
sample  system  studied. 
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1.  INTRODUCTION 


The  seminal  paper  of  Linin'  provides  a  complete  and  definitive  analysis  of  the  structure  of  quasi- 
one-dimensional  diffusion  flames  with  a  one-step  irreversible  reaction.  The  flame  structure  and 
response  are  found  to  be  mainly  described  by  a  Damkohler  number.  Da,  defined  as  the  ratio  of  the 
characteristic  flow  rime  to  the  characteristic  reaction  time.  The  Burke-Schumann  flame  sheet  is 
attained  in  a  distinguished  limit  of  infinite  Da.  With  continuous  decrease  in  Da,  the  characteristic 
extinction  turning-point  behavior  is  exhibited,  yielding  a  distinct  lower  limit  in  Da  beyond  which 
steady  burning  is  not  possible.  This  limit  is  therefore  identified  as  the  extinction  Damkohler  number, 
DaE  K.  The  cause  of  extinction  is  excessive  reactant  leakage  through  the  flame. 

The  system  analyzed  by  Linan  is  a  conservative  one  in  that  in  the  flame-sheet  limit  the  flame 

* 

temperature  is  the  adiabatic  flame  temperature,  Tf.  The  analysis  therefore  neglects  an  inherent 
nonadiabatic  feature  of  flame  processes,  namely  radiative  heat  loss  from  the  flame.  While 
quantitatively  it  is  clear  that  radiative  heat  loss  is  expected  to  reduce  the  flame  temperature  and 
promote  extinction,  there  is  also  a  subtle,  qualitative  feature  which  has  not  been  adequately  recognized 
and  studied.  That  is,  since  radiative  heat  loss  is  a  volume  phenomenon,  the  extent  of  loss  and  thereby 
reduction  in  the  flame  temperature  are  expected  to  increase  with  increasing  system  dimension.  Since  Da 
frequently  also  increases  with  increasing  system  dimension  and  thereby  residence  time,  the  effect  of 
radiative  heat  loss  should  then  increase  with  increasing  Da.  When  such  a  loss  becomes  excessive,  the 
reduction  in  flame  temperature  could  again  cause  extinction.  Thus  it  is  possible  that,  in  addition  to  the 
existence  of  a  lower  limit  in  Da,  namely  DaE  K,  which  characterizes  extinction  of  a  nonradiative  or 
slightly  radiative  flame,  there  could  also  exist  an  upper  limit  in  Da,  say  DaE  R,  which  characterizes 
the  maximum  amount  of  heat  loss  the  flame  can  sustain  for  steady  burning,  and  therefore  can  be  termed 
as  a  radiative  extinction  Damkohler  number.  The  subtlety  here  is  that  since  chemical  reactivity  and 
radiative  loss  both  increase  with  increasing  Da,  it  is  not  clear  a  priori  whether  and  when  radiative 
loss  can  be  sufficiently  strong  to  dominate  over  chemical  enhancement  to  cause  extinction. 

j 

The  possibility  that  an  upper  and  a  lower  bound  in  Da  can  exist  was  suggested  by  Tien  and 

•a 

Sibulkin  from  their  numerical  solutions  of  the  stagnation  and  free-convection  diffusion  flames 
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respectively.  Radiative  loss  in  these  cases  is  from  the  surface.  The  problem  of  diffusion  flame 
extinction  with  flame  radiation  was  analyzed  by  Sohrab  et  al.4  for  the  counterflow  configuration  bv 
using  activation  energy  asymptotics.  A  criterion  for  flame  extinction  was  presented,  although  the 
extinction  turning  point  behavior  of  interest  here  was  not  explicitly  displayed.  Furthermore,  while 
reactant  leakage  is  essential  for  diffusion  flame  extinction,  with  or  without  radiative  heat  loss,  this 
analysis  has  suppressed  reactant  leakage  by  assuming  complete  reaction. 

The  objective  of  the  present  study  is  therefore  to  provide  a  rigorous  analysis  of  the  structure  and 
extinction  of  diffusion  flames  with  flame  radiation.  The  results  demonstrate  conclusively  the  existence 
of  the  dual  extinction  boundaries  of  DaE  R  and  DaE  R.  The  analysis  is  conducted  via  the  model  problem 
of  droplet  combustion  for  its  simple  configuration  as  well  as  ready  conceptual  identification  of  either 
the  droplet  or  the  flame  diameter  as  the  characteristic  system  dimension. 

2.  FORMULATION 

2.1  Governing  Equations 

The  problem  to  be  analyzed,  shown  as  the  inset  in  Fig.  1,  is  the  quasi-steady,  isobaric, 
spherically-svmmetric  droplet  combustion,  with  second-order  Arrhenius  kinetics,  in  a  given  oxidizing 
environment.  Overall  mass  conservation  yields  the  constant  mass  evaporation  rate,  m  =  47tr* pu,  where 
r  is  the  radial  distance,  p  the  density,  and  u  the  radial  velocity. 

For  the  radiative  heat  loss,  recognizing  that  the  heat  loss  rate  per  unit  volume,  qR,  is  a  highly 
sensitive  function  of  temperature,  following  Sohrab  et  al.4  we  approximate  qR  by  an  Arrhenius-type 
function 

qR  =  4a>cT4*BR  exp{  -  Ta  R  /  T )  , 

where  cr=1.36xl0'8  cal/m2-sec-K4  is  the  Stefan-Boltzmann  constant,  k  the  Planck’s  mean  absorption 
coefficient,  BR  a  pre-exponential  rate  parameter  and  Ta  R  an  "activation  temperature"  for  radiative 
loss. 

With  conventional  transport  property  assumptions  (see,  for  example.  Ref.  5),  the  nondimensional 
conservation  equations  for  energy  and  species  are 
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1  d 


_•>  df  ) 


—  df  j  =  Da(p/TfrYFY0  exp(-TaK/T)  -  Ra  exp(  -  Ta  R  /  T  )  , 


dr 


.7 


lhYi_r_:j=_Da(p/TfPYFYoexp(-TaX/T)  -  i  =  F.O  , 
with  the  boundary  conditions 

r=l  T=TS  ,  dT  /  dr=  rhqv  ,  m  Yp  -  (dVF/dr)  =  m  ,  m  Y0  -  (dY0/dr)  =  0 
r 00  :  ,  Yf  -»0  ,  Y0  -»Y0 <00  . 

The  vanous  nondimensional  quantities  are  given  by 


T  = 


T 


qF^ 


'  yf  -  Yp  -  Yq  -  Y0  / 


vO^O 
VFVVF  J 


,  r=r/rs  ,  m  = 


m 


4  ji  rc  X 


/cP 


/  P  =  p/Pf 


qv=qv/qF  -  Ta,K: 


1  a,K 

qF^ 


,  T 


T 


a,R: 


a,R  ^  vO  Bk  P2  cp3  rs2  n  BRrs2S 
,  Da  - r-4 —  -  Ra  = - - 


qF/Cp 


WFqF2R2X 


XqF 


(2.1) 


(2.2) 


(2.3a ) 
(2.3b) 


where  the  subscripts  s  and  f  respectively  denote  quantities  at  the  droplet  surface  and  flame  sheet.  In 
the  above  T  is  the  temperature,  Y;  the  mass  fraction  of  species  i,  X  the  thermal  conductivity,  cp  the 
specific  heat  at  constant  pressure,  D  the  mass  diffusion  coefficient,  qF  the  heat  of  combustion  per  unit 
mass  of  fuel,  v,  the  stoichiometric  coefficient,  W,  the  molecular  weight,  BK  the  pre-exponential  factor, 
Tj.k  th?  activation  temperature,  p  the  pressure,  R  the  gas  constant,  qv  the  latent  heat  of  vaporization 
at  temperature  Ts,  and  the  subscripts  F  and  O  respectively  denote  fuel  and  oxidizer.  The  equation  of 
state  pT=p/R  is  also  applied. 

Using  Eq.  (2.2),  a  coupling  function  (YF-Y0)  can  be  obtained,  yielding  an  explicit  relation 
between  YF  and  Y0 , 


Y0  =  YF -l  +  (  1 +Y0oo)exp(-m/r)  .  (2.4) 

The  problem  is  now  reduced  to  solving  the  mass  burning  rate  m  as  well  as  the  remaining  conservation 
equations  for  YFand  T. 

2.2  Discussions  on  Limiting  Situations 

Several  limits  can  be  discussed  based  on  the  values  of  Ta  K,  Ta  R,  Da  and  Ra.  First,  in  the  limit  of 
f*a,R ~ w'th  Ra  fixed,  we  retrieve  the  adiabatic  droplet  combustion  case.  Similarly,  for  Ta  R— and 
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Da  fixed,  we  have  pure  droplet  vaporization  in  a  radiative  environment.  The  simultaneous 
requirement  of  both  of  the  above  limits  yields  the  case  of  droplet  vaporization  in  a  nonradiative 
environment. 

We  next  take  the  distinguished  limit  of  Ta  and  Da-»°°  for  the  reaction  zone.  For  Ta  R— 
with  Ra  fixed,  radiative  heat  loss  is  frozen  and  we  retrieve  the  classical  adiabatic  droplet  flame- 
sheet  combustion  limit.  If  we  in  addition  take  the  distinguished  limit  Ta  R  — ><»  and  Ra— for  the 
radiative  heat  loss,  then  we  obtain  the  nonadiabatic  flame-sheet  situation.  The  problem,  however,  is 
not  closed  at  the  limit,  as  demonstrated  in  the  Appendix.  Thus  broadening  of  the  radiation  zone  is  an 
inherent  feature  of  a  radiative  diffusion  flame. 

We  next  consider  situations  in  which  one  or  both  of  Ta  K  and  Ta  R  are  finite.  Since  both  the 
reaction  and  radiative  loss  processes  are  temperature  sensitive,  we  have  Ta  K  »  1  and  TaR»l. 
Concerning  the  relative  orders  of  Ta  K  versus  Ta  R,  the  situation  of  Ta  R»  Ta  K  is  of  no  interest  because 
radiative  loss  is  frozen.  On  the  other  hand  the  reverse  situation  of  Ta  K  »  Ta  R  » 1  is  clearly  of  physical 
relevance.  Here  a  thin  reaction  zone  of  0(0  thickness  is  embedded  within  a  somewhat  thicker,  0(6), 
radiation  zone  which,  however,  is  still  much  thinner  than  the  entire  flow'  field  of  dimension  0(1),  as 
shown  in  the  inset  of  Fig.  1.  This  is  the  situation  to  be  analyzed  in  the  next  section. 

The  situation  of  Ta  R  =  0(Ta  K)  is  not  studied  because  it  is  analytically  more  complex,  as  pointed 
out  in  Sohrab  et  a!.4  Its  behavior  is  expected  to  be  intermediate  between  those  of  the  Ta  K»  Ta  R  and 
Ta  K«  Ta  R  situations,  representing  the  transition  between  them.  Details  of  the  transition,  however, 
could  be  interesting  and  merit  further  study. 

We  now  present  the  asymptotic  analysis  of  droplet  combustion  with  flame  radiation. 

3.  ASYMPTOTIC  SOLUTION 

3.1  Outer  Solution 

In  the  limit  of  Ta  K  »  1,  chemical  reaction  occurs  only  in  a  thin  region  of  Ofe)  around  the  location 
of  maximum  temperature.  In  the  outer  regions  bounded  away  from  this  reaction  zone,  chemical  reaction 
is  frozen  because  of  the  reduced  temperature.  Since  the  flow  here  is  influenced  by  two  singularities  of 
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Oie)  and  0(5)  respectively,  the  properties  must  be  expanded  through  both  e  and  5.  Thus  the  proper 
expansion  of  the  properties,  for  example  the  mass  burning  rate  m,  is 

•y  y  "> 

m-  [  m^  ~  e  rh|  +  Ot  e"  )  1  +  5  i  m2  +  e  m3  +  CX  ) )  +  0(  5  ‘  )  ,  (3.1) 

where  the  0(em)  and  0(5n)  terms  are  the  corrections  due  to  finite-rate  kinetics  and  radiative  heat  loss 


respectively. 

In  the  presence  of  radiative  heat  loss,  the  outer  solutions  of  Ypand  T  do  not  take  the  same  form 
and  therefore  need  to  be  separately  analyzed. 

3.1a  Solution  for  Yr  -  Substituting  the  expansion  of  Yp,  in  a  form  similar  to  that  of  Eq.  (3.1),  into 
thi  nonreactive  limit  of  Eq.  (2.2),  and  solving  subject  to  the  boundary  conditions  in  Eq.  (2.3),  we  obtain 


1  <  r  <  f(  :  Y p  q  =  1  -  by  exp(  -  my  /  r )  , 

Yp~  =  - 1  b’  -  |  by  m,  /  ?  )  1  exp(  -  my  /  r )  ,  j  =  l,2; 

rf  <  r  <  °°  :  Ypy  =  by  [  1  -  exp(  -  rhy  /  r )  1  , 

Yf+(  =  bjf-  [  b^  -  ( bymj/r )  ]  exp(  -  my  /  r)  ,  j  =  1 , 2  , 


(3.2a) 


(3.2b) 


(3.2c) 
(3. 2d) 


where  the  b's  are  the  integration  constants  to  be  determined. 

It  is  important  to  recognize  that  while  heat  loss  is  not  explicitly  present  in  the  conservation 
equation  for  Yp,  Yp  is  affected  by  heat  loss  through  its  influence  on  the  mass  burning  rate  rh. 

3.1b  Solution  for  T  -  In  order  to  properly  solve  for  the  temperature  distribution  in  the  outer 
regions,  we  need  to  first  separately  analyze  the  radiation  and  radiatively-frozen  zones  and  then  match 
their  solutions.  In  the  frozen  zone,  energy  transport  is  characterized  by  convection-diffusion  balance. 
The  temperature  fields  are  expanded  in  the  manner  of  Eq.  (31)  which,  when  substituted  into  Eq.  (2.1) 
and  solved  subject  to  the  relevant  boundary  conditions  in  Eq.  (2.3),  yields 


l<r<rf:  t^  =  Ts-qv  +  qvexp{  myf  1  - ( 1  /  r)  ]  }  ,  (3.3a) 

T7=qvm)[l-(l/r)]exp{m0[l-(l/f)]}  ,  j  =  l,2  ;  (3.3b) 

rf<r<°°:  T$  =  T*,  +  ay  ( 1  -exp( -my/r)  ]  ,  (3.3c) 

T|  =  a[-[a[-(aomj/r)]exp(-m0/r)  ,  j  =  l,2,  (3.3d) 

where  the  a's  are  the  integration  constants  to  be  determined. 
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In  the  radiation  region  of  0(5),  the  coordinate  is  stretched  as  C=(r-7f)/8  while  the  proper 
temperature  expansion  is 


T  -  Tf  -  5  [  0]  -  ( e  /  5  )  02  -  e  ©3 


■a  5“  )  . 


Substituting  Eq.  (3.4)  and  £  into  the  cherrucallv-frozen  form  of  Eq.  (2.1)  and  expanding,  we  obtain 


d'©f  /ds  =  -  Ra  exp(  -  ©7  )  , 

d'©“  / d^~  =  -  Ra  ©r  exfx  -  ©j  )  ,  j  =  2 , 3  , 


(3.5a) 

(3.5b) 


where  the  small  expansion  parameter  5  is  defined  as  8=Tf /Ta  R  while 

Ra  =  8  Ra  expt  -  Ta  R  / T. )  (3.6) 

is  a  reduced  parameter  representing  radiative  heat  loss.  Thus  energy  conservation  in  the  radiation 
region  is  diffusion-radiation  balanced. 

Matching  the  solutions  in  the  radiation  and  frozen  regions  as  ?-*  rf  yields 

Tf  =  Ts-qv  +  qvexp{m0[l-(l  /  rr )  ] }  ,  a£  =  ( Tf-T*,)  /  [  1  -exp< -mp/ r( )  ]  ,  (3.7) 

and  the  boundary  conditions  needed  to  solve  Eqs.  (3.5)  as  £-++<»,  giving 


C<0  : 

©1'  =  21nil  -  [Ra/(2g5)]exp(g0i;  +g,  )l  -  gpC-g,  , 

€6  =  -  ( Tf  -  Ts  +  qv )  ( ih)  / go )  [  1  -  ( 1  /  r) )  ]  [  go  +  2  Ra  exp(  -  SJ  )  ]1/2 


1  1  3  1  +  (Ra/g0)Cexp(gbC  +  gi)-{[Ra/(2go)]exp(g0C  +  g1)} 


— —  =  -gnm,  1  +  — - — 

d^  60  'l  r, 


{1  -  [Ra /(2g^)]exp(g0i;  +  g1)}2 


C>0  : 

©^  =  2  In  1 1  -  [  Ra  /(  2  h^  )  ]  exp(  h,  -  ho  £  )  I  +h0  C  -  hj  , 

©£=  -{ (  aj  /  ho )  1 1  -  exp(  -  m0/ff)]  +  (m1ff/m0)}[h^  +  2Ra  exp(  -  ©j )  ]I/2 

d©3  .  2 

-^r-  =  (a^mo/rj  )exp<-m0/rf  i  +  hpmj  [(1  /  mo)-(l  /  r£ )  1  }x 

1  -  (  Ra / hp )  ^  exp(  h,  - h0 Q -{ [  Ra /(  2  hp  )  ]  exp(  hj-hpC)}2 

{l-[Ra/(2h20)]exp(h1-h0O}2 


(3.8a) 


(3.8b) 


(3.8c) 


(3.9a) 

(3.9b) 


(3.9c) 


where 
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g0  =  (Tf-Ts-qv)m<)/rf‘  ,  l^  =  ( Tf  -  )  (  m^,  /  rf2  )  exp(  -  mg/  rf )  /  [  1  -  exp<  -  mg  /  r( )  ]  , 

g1=(Tf-Ts~qv)m;[l-(l/r{)]  ,  h, -a,  (  1 -cxp( -mo/Pf )  J  +  (  h0  ff  rh2/if^, )  .  (3.10) 

3—  Inner  Solution 

In  the  reaction  region  of  O(e)  thickness,  the  temperature  is  within  an  O(e)  reduction  from  the 
maximum  temperature  Tf.  The  stretched  spatial  coordinate  is  ^  =  (7 -rf)/e  while  the  temperature  and 
fuel  distributions  are  expanded  as 


f  =  i,-e[  0j+O(  5)*O(c/5)-t----]-e2102  +  o(5)  +  O(e/S)  +  -"]  +  o(  e3 ) 
Yp~£[<f>]  +  Ot  5  )  -^  Q  £  /  5  )  +■■■!  +  [  <>2  +  Of  5  )  +  Of  £  /  8  )  +  •■■]  +  Of  £  )  . 


(3.12) 


Since  oxidizer  concentration  is  also  an  0(£)  quantity  in  this  thin  region,  substitution  of  Hq.  (3.12)  into 
Eq.  (2.4)  and  forcing  the  leading  order  term  to  vanish  yields 


rf  =  rh0  /  lnt  1  +  Y0  x  )  . 


(3.13) 


Substituting  Eqs.  (3.11)  and  (3.12)  as  well  as  the  stretched  coordinate  \  into  Eqs.  (2.1)  and  (2.2), 
and  expanding,  we  have  two  source-free  equations  expressing  the  local  coupling  functions  and  a  third 
equation  with  the  source  term  representing  chemical  reaction. 


d2 

—2  (  ©I  -  0i  )  =  0 

a; 


2  d  T  d 


- -  (  -  db  )  +  t - pr  ,  £,  —  (  0j  —  )  —  ~Tr  (  9]  -  0; 

d\-  rt  d^~  J  V 


d 


-0i)  =  0  , 


(3.14) 


(3.15) 


-T  =  Da  p,  [<>,-(  m,  /?()  +  (  riio  /  rf‘  )  ^  1  exp(  -  8,  )  , 


(3.16) 


where  the  small  parameter  e  is  defined  as  c  =  Tf /Ta  K  and 
Da  =  e3  (  Da  /  T, 2  )  exp(  -  Ta  K  /  Tf ) 


(3.17) 


is  a  reduced  Damkohler  number.  In  obtaining  Eq.  (3.16),  we  have  also  used  p  =  1  +0(e). 

Equations  (3.14)  to  (3.16)  show  that  this  region  is  reactive-diffusive.  Radiative  heat  loss  is  not 
important  here  because  its  volume  is  much  smaller  than  the  total  volume  of  the  radiation  region. 

33  Final  Solution 


The  boundary  conditions  required  to  solve  Eqs.  (3.14)  to  (3.16)  are  derived  from  matching  the  inner 
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and  outer  solutions  as  r  -»rf  for  YFand  £->0  for  T.  Upon  matching  we  first  obtain  bg  =  1  +Y0oo/  bQ=0  and 
the  solutions  for  rii,,  aj and  b?.  Next,  the  boundary  conditions  for  the  inner  equations  as  are  also 

obtained.  Integrating  Eq.  (3.14)  twice  and  Eq.  (3.15)  once,  and  applying  the  proper  boundary  conditions, 
we  obtain 

(  go  2  Ra  )  ^  +  (  +  2  Ra  )  =  m^  /  r^  ,  (3.18) 

and  two  additional  relations  to  express  b^  in  terms  of  a|  and  fn,.  With  these  relations,  Oj  can  be 
expressed  in  terms  of  0]  as 

di  =  {  0,  +  [rft,  /  ln(  1  +  Y0  oo)  1  }  (  h^  +  2  Ra  )1/2  .  (3.19) 

Substituting  Eq.  (3.19)  into  Eq.  (3.16)  yields  a  second  order,  nonlinear  ordinary  differential  equation 
with  0|  as  the  only  unknown,  subject  to  the  matching  conditions  obtained.  This  equation  needs  to  be 
solved  numerically. 

The  problem  can  h_  _•  oidered  to  be  completely  solved  at  this  stage  in  that  all  parameters  to  the 
lowest  orders  needed  io  solve  the  problem  are  either  determined  or  are  determinable  through 
straightforward  numerical  integration.  The  solution  includes  three  expressions,  Eqs.  (3.7),  (3.13)  and 
(3.18),  whvh  define  the  three  bulk  combustion  properties  of  the  flame  temperature  Tf,  the  flame-sheet 
location  rf,  and  the  mass  burning  rate  ir^,  in  terms  of  the  system,  kinetic  and  radiation  parameters. 
Furthermore,  solution  of  the  differential  equation  Eq.  (3.16),  with  (Jq  given  by  Eq.  (3.19),  determine  the 
parameters  a|  and  rh,  and  consequently  the  condition  for  flame  extinction. 

4.  RESULTS  AND  DISCUSSIONS 

4.1  Re-Scaling 

The  solutions  just  obtained  are  not  amenable  for  ready  interpretation  of  the  system  behavior 
because  the  quantifying  parameters  8,  e,  Ra  and  Da  are  all  functions  of  the  leading-order  flame 
temperature  Tf,  which  in  turn  depends  on  the  characteristics  of  radiative  heat  loss.  It  is  therefore 
necessary  to  re-scale  all  the  parameters  to  a  fixed,  absolute  reference  state.  Failure  to  do  so  could  lead 
to  spurious  interpretation  of  the  results. 

The  obvious  choice  of  the  reference  state  is  the  adiabatic  state.  Designating  properties  of  this 
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state  by  the  superscript  *,  we  have 


T*  =[Tx~  Y0,oo(  1  -qv-t,)] /(  1  -  YQso)  ,  (4.1) 

hTo=  ln{  1  -  [  (  Y0oo-Tb )  /qv]  }  ,  (4.2) 

rf  =11^/  ln(  1  +  Y0  =o  )  .  (4.3) 

The  parameters  referred  to  this  adiabatic,  reference  state  can  then  be  expressed  as 
Da  =  e  3  Da  exp(  -  Ta  K  /Tf  )  /Tf  “  ,  e  =Tf  ~/Ta  K  , 

Ra  -  8  Ra  exp(  -  Ta  R/Tf  )  ,  8=Tf2/TaR  .  (4.4) 

With  the  expressions  in  Eq.  (4.4),  the  system  dependent  parameters  can  be  rewritten  as 

Da=  Da’(  %/j]  )4exp[  Ta>  K  (  t|  _1  - ff1 )  ]  ,  (4.5) 

Ra  =  Ra‘(  Tt/Tf  )2exp[Ta  R(Tf  “’-tf'1)]  ,  (4.6) 


■f  ~  *  ~  ")  +  r— 

while  aj  =(Tf/Tf  )'a,  is  obtained  from  ea[  =e  a:  and  a  similar  expression  can  be  obtained  for  mj . 

—  *  ~  # 

Since  Da  and  Ra,  and  hence  Da  and  Ra  ,  vary  with  the  droplet  size  in  the  same,  quadratic, 
manner,  we  can  define  a  parameter  F  =  Ra  /Da  ,  which  is  the  ratio  representing  the  characteristic 
reaction  time  to  the  characteristic  radiation  time.  It  depends  only  on  the  kinetic  and  radiation 
constants  of  the  system. 

In  addition  to  the  need  for  re-scaling  to  an  absolute  reference  state,  it  is  also  necessary  to  clearly 

identify  the  roles  of  the  various  dependent  and  independent  parameters  in  their  dimensional  and 

nondimensional  forms.  Of  particular  importance  is  the  droplet  size,  r^  which  has  been  extensively  used 

in  nondimensionalization.  For  example,  an  increase  in  the  (reduced)  Damkohler  number  Da  can  be 

interpreted  as  either  an  increase  in  chemical  reactivity  or  an  increase  in  the  droplet  size.  Similarly,  an 

—  » 

increase  in  the  heat  loss  parameter  Ra  can  imply  either  an  increase  in  the  loss  intensity  or  an  increase 

in  the  droplet  size.  Thus  a  clear  measure  of  the  relative  influence  of  reaction  and  radiation  is  the 

parameter  T,  in  which  the  droplet  size  is  divided  out.  The  problem  is  defined  by  using  any  two  of  the 
—  *  —  * 

three  parameters  Da  ,  Ra  and  T. 

It  is  further  found  that  the  final  structure  equation  of  the  present  system  can  be  converted  to  that 
of  Lifian1  by  re-scaling  the  various  parameters  and  variables  as 
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0  =  aI/3  0[  -  yq  -  a  =  4  Da*  (  Tf  /Tf  )4  [rf  /  ln<  1  +  Y0oo)  ]  2exp[Ta  K  (  Tf  1  - Tf  1  )  ]  , 

T1  =  (  a1  n  /  2  )  i  4  !n(  1  +  \Qx  i  -  m,  |  /  r,  ,  y  =  1  -  l  2  rf  /  ln(  1  +  Y0oo  )  ]  (  g^  ^  2  Ra  )W2  #  . 
The  structure  equation,  Eq.  (3.16),  is  then  reduced  to  that  of  Linan,  given  by 


d20  -2  •>  _i/3- 

—  =  (92-Ti‘)exp[-a  1/3  (  0  +  yn 
dTf 


dO  ) 


d0 

dff 


=  1 


(4.7) 


The  quantities  a  and  m,  can  be  determined  from  solving  Eq.  (4.7)  as 


a  =  -  a^  *  Y0oo  /  (  1  +  Y0oo  )  =  a~  1  73  ( Tf  /  Tf*  )2  [  /  (  h^  +  2  Ra  /  22  ]  lim  (  0  - 1| )  ,  (4.8) 

T]  -40 O 

m,  =- {  a'  1  73  In  (  1  +  Y03O  )  /  [  rt  (  g2- 2  Ra  )’  12  )  }  (  Tf /Tf  y  lim  (0  +  q)  .  (4.9) 

T]  — >  -  oo 

Equation  (4.7)  has  been  numerically  solved  by  Linan.1  Results  show  that  extinction  is  expected  for 
u<ac,  where  an  is  approximately  given  by 

aE  =  e  [  ( 1  -  ! yl  >  -  (  1  -  !y!)2  +  0.26  (1-  ly!)3  +  0.055  (  1  -  ! y I  )4  ]  (4.10) 

with  very  good  accuracy.  It  may  be  noted  that  in  the  presence  of  radiative  heat  loss,  y  varies  with  a 
through  Da  for  fixed  T  and  therefore  cannot  bo  independently  specified. 

We  mention  in  passing  that  Sohrab  et  al.4  assumed  complete  reaction  so  that  there  is  no  reactant 
leakage  through  the  flame.  This  means  that  in  addition  to  the  required  boundary  conditions  to  solve 
Eq.  (4.8),  the  additional  constraints  (0 +  q)_oo=(0-T))oo  =  O  are  applied,  which  renders  the  problem  over¬ 
determined. 

4.2  Results 

To  demonstrate  the  salient  features  of  the  solution,  the  combustion  of  a  heptane  droplet  at  its 
boiling  temperature  and  burning  in  the  standard  atmosphere  is  solved.  The  system  parameters  adopted 
are 


Ts  =  372  K  ,  Too  =  298  K  ,  Ta  k  =  22,000K  ,  TaR  =  8,000K  ,  Cp  =  0.317  cal /gm-K  , 
qp  =  10.8 kcal  / gm  ,  qv  =  75.8 cal  / gm  ,  (v0W0)/ (vFWF)  =  352  ,  Yo,c»  =  0232  . 

We  first  present  the  flame  characteristics  when  the  reaction  zone  is  infinitely  thin  and  reactant 
leakage  is  suppressed.  The  flame  is  still  of  finite  thickness  due  to  radiation  broadening.  The 
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combustion  parameters  T(,  ri^,  and  rf  are  given  by  the  three  algebraic  expressions  Eqs.  (3.7),  (3.13)  and 
'3. IS),  which  can  be  iteratively  solved  together  with  the  scaling  relation  Eq.  (4.6). 

Figure  1  shows  the  flame  temperature  T(  and  mass  burning  rate  n\,  as  functions  of  Ra*.  The  flame 

radius  7,  is  not  shown  because  it  is  directly  proportional  to  itiq.  It  is  seen  that  T(  and  rry,  both  decrease 

~  * 

with  increasing  Ra  .  Thus  for  a  fixed  droplet  radius  rs,  increasing  radiative  loss  intensity  reduces  the 
flame  temperature  and  mass  burning  rate,  as  expected.  However,  if  we  interpret  the  increase  in  Ra"  as 
due  to  an  increase  in  rs,  with  fixed  radiation  parameters,  then  we  should  also  account  for  the  fact  that 
the  physical  burning  rate  rr^  varies  linearly  with  rs.  Then  the  numerical  values  show  that  my,  still 
increases  with  rs,  albeit  at  a  slower  rate  in  the  presence  of  radiative  loss. 

Since  the  reaction  rate  is  assumed  to  be  infinitely  fast,  extinction  is  not  described  at  this  level  of 
solution.  Consequently  Fig.  1  does  not  exhibit  the  extinction  turning  point  behavior.  To  allow  for  finite- 
rate  kinetics,  Eq.  (4.7)  is  solved  numerically.  Figure  2  plots  a  versus  Da  for  several  values  of  Ra  , 
where  a  represents  the  temperature  reduction  due  to  incomplete  reaction  and  is  directly  related  to  the 
amount  of  fuel  leakage  through  the  flame. 

We  first  discuss  the  adiabatic  limit,  Ra  =0,  which  corresponds  to  y=  0.867  for  the  present  system. 
This  is  the  situation  studied  by  Linan.1  It  is  seen  from  Fig.  2  that  there  exists  a  minimum  Da  ,  Da^, 
such  that  there  is  no  solution  for  Da  <  DaEK  and  there  are  two  solutions  for  Da'  >  Da^K.  Thus  Da^K  can 
be  identified  as  a  (reduced)  kinetic  extinction  Damkdhler  number.  The  lower  branch,  which  shows  a 
decreases  with  increasing  Da  ,  is  the  physically  realistic  solution. 

To  investigate  the  influence  of  radiative  heat  loss,  it  is  important  to  specify  the  process/ 
parameter  which  is  held  fixed  in  the  comparison.  To  demonstrate  this  point,  in  Fig.  2  the  effect  of 
radiation  is  exhibited  by  plotting  a  versus  Da  on  constant  Ra  .  Since  both  Da*  and  Ra*  depend  on  the 
droplet  size  through  r9,  rs  should  be  considered  as  a  fixed  quantity  in  the  comparison.  Figure  2  then 
shows  that,  for  a  given  reference  reactivity  Da  ,  increasing  radiative  loss  through  Ra’  increases  the 
extent  of  reactant  leakage.  This  is  reasonable  because  of  the  reduction  in  the  flame  temperature  and 
thereby  reaction  rate  due  to  the  heat  loss.  The  kinetic  extinction  Damkohler  number  is  consequently 
also  increased.  By  the  same  reasoning,  we  see  that  for  a  given  radiative  loss  Ra  ,  increasing  reactivity 
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IX)  reduces  the  extent  of  leakage. 

I tgure  2  is  relev, mt  tor  studying  the  efforts  of  changing  the  fuel /oxidizer  system  <e  g  luel 
substitution,  oxidizer  enrichment)  on  a  droplet  ot  fixed  si/e.  However,  frequently'  a  more  rele,  mt 
cpiestion  to  ask  is  the  effect  of  varying  the’  droplet  size  for  a  fixed  fuel/oxidizer  system.  For  example, 
after  a  fuel  droplet  or  particle  has  achieved  ignition,  the  droplet  size  rs  will  steadily  decrease 
although  its  physical  flame  size  r,  can  first  increase,  due  to  fuel  vapor  accumulation,  ’  and  then 
decrease  The  relevant  plot  to  demonstrate  droplet  size  effects  is  again  a  versus  Da  ,  but  on  constant  F, 
as  shown  in  Fig.  3.  For  such  a  plot  we  can  consider  the  reference  chemical  reactivity  to  be  fixed.  Thus 
lor  a  given  Da  ,  inc  reasing  I  increases  the  influence  of  radiation,  while  for  a  fixed  F  increasing  Da 
implies  increasmg  values  of  rs. 

Figure  3  shows  that,  for  smaller  values  of  Da  ,  the  flame  response  is  qualitatively  similar  to 

that  of  Fig.  2,  exhibiting  the  turning  point  behavior  at  DaL:R.  The  behavior,  however,  becomes 

qualitatively  different  for  larger  values  of  Da  .  It  is  seen  that  by  continuously  increasing  Da  along  the 

lower  branch  of  a  constant  F  curve,  the  reactant  leakage  first  decreases  and  then  increases  until  another 

turning  point  is  encountered.  According  to  our  earlier  discussion  and  anticipation,  this  turning  point 

should  represent  an  alternate  mode  of  extinction,  namely  extinction  due  to  fuel  leakage  in  the  presence 

of  excessive  radiative  heat  loss  and  (lame  temperature  reduction  for  large  droplet  sizes.  We  therefore 

designate  this  critical  state  as  that  of  the  radiation  extinction  Damkbhler  number,  Da^R.  It  may  be 

emphasized  that  this  mode  of  extinction  is  caused  by  the  slow  kinetics  due  to  the  heat  loss  but  not  the 

~  •  'w  •  ~  •  "*'• 

loss  itself.  Thus  steady  burning  is  possible  only  for  Da^K<  Da  <  Da^R.  This  range  in  Da  is  further 
diminished  with  increasing  T  (Fig.  3)  such  that  there  exists  a  maximum  T,  say  rc,  beyond  which  burning 
is  absolutely  not  possible,  for  the  fuel/oxidizer  system,  regardless  of  the  droplet  size. 

In  Fig.  4  we  have  plotted  the  extinction  Damkohler  number  DaE  and  the  corresponding  flame 
temperature  as  a  function  of  F.  This  curve  can  be  generated  directly  by  using  Eq.  (4.10)  with  the  proper 
definitions  of  a  and  y.  Thus  for  a  given  T,  the  possible  values  of  Da  are  bounded  by  the  upper  and 
lower  branches  of  this  curve,  which  respectively  represent  the  kinetic  and  radiative  extinction 
Damkohler  numbers  Da^K  and  Da^R.  The  maximum  value  of  T  is  found  to  be  fc=0.0302  for  this  model 
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system,  it  is  also  seen  that  the  flame  temperature  reduction  is  small  for  the  upper,  kinetic  extinction 
regime,  but  is  substantially  larger  for  the  lower,  radiative  extinction  regime,  as  is  reasonable  to  expect. 

It  is  of  interest  to  estimate  the  characteristic  radiative  extinction  droplet  size.  Such  an 
estimation  requires  values  of  a.  and  BR.  For  a.  we  used  2.97xl0'4  cal/cm-sec-K,  which  is  that  of  air  at 
2,000K.  To  estimate  BR,  we  first  fitted  T4  in  the  temperature  range  of  1,600K<T<2,500K  to  yield 
T4  =  8.90xlOuexp(-8,OOOK/T)  K4.  Next,  the  Planck's  mean  absorption  coefficient  is  taken  as  the  molar 
averaged  value  of  the  mixture  in  the  radiation  region.  Since  the  amount  of  reactants  is  small  compared 
to  that  of  the  products  in  this  region,  it  is  reasonable  to  consider  that  the  mixture  consists  of  C02,  H20 
and  NS .  Furthermore,  since  \N  is  very  non-radiative  as  compared  to  C02  and  H-,0,  we  have 
k  =  xco2Kco2+xH20KH20-  Using6  kCo2  =  5.63x10'2  cm'1  and  kH2q  =  7.13x10‘3  cm’1  evaluated  at  1  atm  and 
2.000K,  and  xco^=0.124  and  xj  ^0=0.142  for  the  stoichiometric  burning  oc  heptane  in  air,  we  obtain 
k  =  8.00x  10'3  cm  1 .  Using  k,  a  and  the  fitted  T4,  BR  is  determined  as  3.87xl03  cal/m3-sec  and 
consequently  Ra  =(28.2  m’2)r2  with  Tf  calculated  from  Eq.  (4.1). 

Experimentally,  n-alkanc  droplets  burning  in  air  typically  extinguish,  in  the  kinetic  limit,  in 

the  range  around  rsEK  =  0.1  mm.  This  yields  Ra^  =  2.82x10'  .  Using  this  value  and  the  relation 

•  —  • 

Da£K=  RaEK/T,  Fig.  4  shows  that  there  is  no  solution  even  with  the  maximum  T  of  0.0302.  iherefore 
the  flame  cannot  exist  at  all.  This  result  obviously  contradicts  experimental  observations  because 
n-alkane  droplets  can  be  routinely  burned  in  room  air.  The  cause  of  this  discrepancy  is  believed  to  be 
the  excessively  large  flame  size  predicted  by  the  theory.  That  is,  it  is  well  established5  that  because 
of  the  unity  Lewis  number  assumption  and  the  neglect  of  the  fuel  vapor  accumulation  process,  the 
classical  d  -law  solution,  which  is  the  leading  order  solution  in  the  present  problem,  predicts 
flamefront  standoff  ratios  for  n-alkane  droplet  combustion  in  air  to  be  typically  40.  Experimentally5 
the  ratio  is  more  in  the  range  of  5  to  10.  This  large  theoretical  flame  size  then  leads  to  excessive 
radiative  heat  loss  and  thereby  promotes  radiation-induced  extinction.  Since  radiative  heat  loss 
scales  with  rf ,  this  effect  can  be  approximately  accounted  for  by  reducing  Ra^K  by  the  square  of  the 
ratio  of  the  experimental  to  the  theoretical  flame  size.  If  we  take  this  ratio  to  be,  say  6,  then  through 
iteration  we  find  Da^K  =  4.67x10'  and  f  =  0.0167,  which  yields  a  radiative  extinction  droplet  size  of 
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rs  E  R=  0.93  mm.  This  value  is  about  10  times  larger  than  the  kinetic  extinction  droplet  size.  It  is 
reasonable  to  expect  that  in  the  presence  of  soot  the  value  of  rs  E  R  could  be  further  reduces  to  render 
radiative  extinction  an  important  process  in  heterogeneous  combustion. 

5.  CONCLUDING  REMARKS 

In  the  present  study  we  have  successfully  analyzed  the  structure  and  extinction  characteristics  of 
diffusion  flames  suffering  radiative  heat  loss  from  the  flame,  when  the  thickness  of  the  radiation  zone 
is  asymptotically  intermediate  of  those  of  the  reaction  zone  and  the  bulk  flow  field.  Perhaps  the  most 
interesting  result  from  the  present  study  is  the  identification  of  the  existence  of  dual  turning  points. 
Specifically,  by  fixing  the  thermochemical  aspects  of  a  fuel/oxidizer  system,  we  have  shown  that 
steady  droplet  burning  is  possible  only  for  droplet  sizes  within  a  certain  range.  Extinction  of  smaller 
droplets  is  caused  by  purely  finite-rate  kinetics,  while  extinction  of  larger  droplets  is  caused  by 
excessive  radiative  heat  loss  which  eventually  again  leads  to  the  kinetic  limitation.  The  range  of 
droplet  size  for  steady  burning  narrows  with  increasing  radiation  intensity,  as  represented  by  the 
parameter  T,  such  that  droplet  burning  is  not  possible  beyond  a  certain  critical  Tc. 

Extension  of  the  present  results  to  other  quasi-one-dimensional  flow  configurations  is 
straightforward,  involving  primarily  changes  in  the  definitions  for  Da  and  Ra  to  account  for  the 
different  expressions  of  residence  times. 
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APPENDIX 

Broadening  of  Radiative  Diffusion  Flames 

With  the  flame  sheet  located  at  rf  and  having  a  temperature  Tf,  the  temperature  and  fuel 

distributions  can  be  solved.  Since  both  the  fuel  and  oxidizer  have  to  be  totally  consumed  at  the  flame 
sheet,  we  obtain  from  Eq.  (2.4)  rf=rh/ln(l  +  Y0 <0o).  Evaluating  the  temperature  at  the  flame  sheet 

yields 

=  "Ts - qv qv cxPi  m  1 1  - ( 1  /  rf )  ]  }  .  (A.l) 

The  jump  relations  can  be  obtained  by  integrating  Eqs.  (2.1)  and  (2.2)  once  across  the  flame  sheet, 
yielding 

?f2[(dtVdr)-  (dT"/cfr)]  =  +  I2  ,  ?f2  (dYp/cfr) -  =- 1,  ,  (A. 2) 

where  the  superscripts  -  and  +  respectively  denote  the  fuel  and  oxidizer  side,  and 


Ij  =  (  Da  /  Tf2  )  f  *  ?2p2  YpY0exp(  -TaK  /  T)  dr  ,  I2  =  Ra  f  '  r2  exp(  -  TaR  / T  ) dr  . 

Jrf 

Substituting  the  solutions  of  T1  and  Yj:  into  Eq.  (2.6)  gives 
m{T«  +  Y0oo-Ts+qv[l-exp(fh)]}  =  Y0oo  I2 


and  Ij=m.  It  is  seen  that  there  are  three  unknowns,  rh,  Tp  and  I2,  which  need  to  be  determined  from  the 
two  relations  of  Eqs.  (A.l)  and  (A.  3).  The  problem  is  therefore  not  closed. 


Figure  1 . 


Figure  2. 


Figure  3. 


Figure  4. 


Variation  of  the  (lame  temperature  T,  and  the  mass  burning  rate  fry,  with  the  radiative  heat 
*'  * 

loss  parameter  Ra  .  The  inset  is  a  schematic  defining  the  problem  analyzed. 

The  extent  of  flame  temperature  reduction  or  reactant  leakage  as  a  function  of  Da  for  fixed 

* 

values  of  Ra  . 

The  extent  of  flame  temperature  reduction  or  reactant  leakage  as  a  function  of  Da  for  fixed 
values  of  f. 

Extinction  Damkohler  number  DaE  and  flame  temperature  Tf  E  as  functions  of  the  radiation 


heat  loss  parameter  f. 
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Abstract 


A  theoretical  analysis  is  performed  to  study  the 
distribution  of  small-size  particles  in  a  flow  adjacent  to  a 
hot  wall.  Because  of  the  thermophoretic  motion,  a  particle 
stagnation  point  is  created  in  the  interior  of  the  flow  which 
can  produce  singularity  when  particle  diffusion  is  neglected. 
By  including  Brownian  diffusion  and  using  matched  asymptotic 
expansions,  the  structure  of  the  inner  diffusion  xayer  near  the 
particle  stagnation  point  is  resolved.  The  presented  results 
illustrate  how  the  diffusion  layer  is  evolved  and  its 
structural  variation  with  a  number  of  parameters. 
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I.  Introduction 

Thermophoresis  is  the  motion  of  small  particles 
relative  to  their  containing  gas  induced  by  a  temperature 
gradient  [1,  2].  Because  of  the  existence  of  thermophoret ic 
velocity,  particle  concentration  distribution  can  have  unusual 
characteristics  which  requires  special  analysis.  In  the 
presence  of  a  cold  wall,  the  thermophoretic  motion  drives  the 
particles  toward  the  wall  which  affects  their  deposition  rate. 
In  most  theoretical  works  concerning  deposition  rates,  e.g., 
Refs.  [3-7],  particle  diffusion  is  neglected.  This  is  because 
the  diffusion  coefficient  is  very  small  which,  in  the  case  of  a 
cold  wall,  will  have  only  negligible  influence  on  the  particle 
convective  flux  toward  the  wall.  The  neglect  of  particle 
diffusion,  however,  lowers  the  order  of  the  differential 
equation  of  the  particle  concentration  equation.  Therefore, 
the  concentration  structure  in  the  thin  layer  next  to  the  wall 
cannot  be  studied.  Walker,  Homsy,  and  Geyling  [8]  included 
Brownian  diffusion  in  their  analysis  of  thermophoretic 
deposition  in  flows  inside  a  cold  tube.  The  structure  equation 
is  solved  using  matched  asymptotic  expansions.  They  found, 
however,  that  because  near  the  wall  the  fluid  velocity  is 
negligible  while  the  thermophoretic  velocity  is  finite,  the 
wall  particle  deposition  rate  is  unaffected  by  the  diffusion 
layer . 
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For  a  hot  wall,  the  thermophoret ic  motion  drives  the 
particles  away  from  the  wall.  When  particle  ui  ^  j.on  is 
neglected,  a  particle  stagnation  point  (or  surface)  is  created 
and  a  dust-free  layer  adjacent  to  the  surface  is  predicted  [3, 
7,  9,  11].  For  low  Mach  number  flow,  Goren  [3]  is  the  first  to 
show  that  when  the  product  of  thermophoretic  coefficient  (K) 
and  Prandtl  number  (Pr)  is  less  than  unity,  the  particle 
concentration  drops  to  zero  quickly  at  the  particle  stagnation 
surface.  But  when  KPr  is  greater  than  unity,  a  singularity 
occurs  -  the  particle  concentration  becomes  unbounded  at  the 
particle  stagnation  surface.  Since  in  both  cases,  the 
concentration  gradients  become  very  steep,  particle  diffusion 
can  become  important  even  when  the  diffusion  coefficient  is 
small.  For  the  case  KPr  <  1,  analysis  [11]  has  been  performed 
for  incompressible  flow  (small  temperature  difference  and  low 
Mach  number)  including  diffusion  to  assess  the  degree  of 
particle  leakage  across  the  particle  stagnation  surface.  But 
for  KPr  >  1,  the  more  severe  case  for  the  solution  misbehavior 
(singularity),  no  remedy  has  been  given. 

In  this  work,  a  theoretical  analysis  is  given  for  the 
particle  concentration  distribution  in  a  low  Mach  number, 
compressible  flow  over  a  hot  wall.  Particle  diffusion  is 
included.  By  using  the  particle  Schmidt  number  as  the  large 
parameter,  an  inner-outer  expansion  technique  is  employed  to 
study  the  particle  diffusion  layer  structure  for  both  the  cases 
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of  KPr  greater  than  unity  and  less  than  unity. 

II.  Governing  Equations 

The  flow  to  be  analyzed  is  the  axisymmetric 
stagnation-point  region  of  a  heated  plate  as  shown 
schematically  in  Fig.  1.  Experimentally,  Phis  configuration 
can  be  created  either  by  impinging  a  jet  perpendicularly  onto  a 
plate  or  by  placing  a  blunt  body  in  a  windtunnel.  The  small 
particles  are  seeded  in  the  upstream.  Fig.  1  shows  the 
existence  of  two  stagnation  points:  a  gas  flow  stagnation 

point  on  the  plate  surface  at  x  =  y  =  0  and  a  particle 

•  • 

stagnation  point  at  x  =  0  and  y  =  y  (v  =  t?  )  .  The  presence  of 
the  particle  stagnation  point  in  the  interior  of  the  flow  field 
is  a  result  of  the  particle  thermophoretic  motion  away  from  the 
heated  plate.  To  analyze  the  structure  of  particle 

concentration  around  this  particle  stagnation  point  is  the 
objective  of  this  work. 


The  equations  governing  the  low  speed  flow  (Mach 
number  square  much  less  than  unity)  in  the  axisymmetric 
stagnation-point  region  are  given  by  the  compressible  boundary 
layer  equations: 


dpux  +  3pvx  m 
9x  dy 


(1) 


=  _  *£  +  L_ 
ax  ay 


(2) 
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(4) 


l  a  ,,  a  „  a 

-  — xpuY  +  — pv  Y  =  — 
x  3x  K  oyK  P  ay 


pD3y 


where  Eqs.  (1)  -  (3)  are  the  conventional  continuity,  momentum 
and  energy  equations  respectively.  The  pressure  gradient  term 
in  Eq .  (2)  is  equal  to  -  p^u^  3u^/5x  and  the  x-component 
free-stream  velocity  is  given  by 


u  =  ax 

e 


(5) 


where  a,  the  velocity  gradient,  is  the  parameter  characterizing 
the  flow.  Eq.  (4)  is  the  mass  conservation  equation  for  the 
particles.  In  Eq.  (4),  Y  is  the  mass  fraction  of  particles,  v 

p 

is  the  particle  velocity  which  is  the  sum  of  gas  velocity  (v) 
and  thermophoretic  velocity  (V  )  : 

v  =  v  +  V  (6) 

p  t  '  ' 

The  expression  for  the  thermophoretic  velocity  is: 

vt  =  -K(M/p)VT/T  (7) 

where  K  is  the  thermophoretic  coefficient  whose  magnitude 
depends  on  the  regime  of  the  flow  [12].  The  value  of  K 
normally  taken  is  between  0.2  and  1.25  and  in  the  present  study 


it  is  treated  as  a  parameter.1  The  last  term  in  Eq.  (4) 
represents  particle  (Brownian)  diffusion  and  D  is  the  diffusion 
coefficient  which  is  normally  very  small. 

The  boundary  conditions  for  the  flow  and  the 
temperature  fields  are: 

y=0:  u=v=0  ,  T=T  (8) 

y-><»:u  =  u,T  =  T  (9) 

and  for  the  particles: 

y  =  0:  Y  =  0  (10) 

y  ~*  °°:  Y  -  Y.  (11) 

The  condition  of  Y  =  0  in  Eq.  (10)  is  the  result  of  assuming 
that  we  have  a  "sticky"  wall  that  will  act  as  a  perfect  sink 
for  the  particles  [15].  The  condition  of  Y  =  Y^  in  Eq.  (11) 
indicates  that  the  particles  are  seeded  in  the  free  stream.  To 
proceed  further,  the  following  assumptions  are  adde^: 

1.  the  gas  obeys  ideal  gas  law  and  the  product  pu  is 
a  constant, 


'Experimental  determination  of  K  is  a  difficult  task  [13].  A 
novel  measurement  technique  using  counterflow  diffusion  flames 
has  recently  been  suggested  [14], 
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2.  gas  heat  capacity  C  ,  the  thermophoretic 
coefficient  K,  and  the  Prandtl  number  Pr  are 
constant, 

3.  the  mass  fraction  of  particles  is  much 
smaller  than  unity  which  does  not  affect  the 
overall  continuity  equation, 

4.  the  Schmidt  number  (based  on  Brownian 
particle  diffusion  coefficient)  is  constant, 

5.  the  particles  are  of  uniform  size  and  inert. 

There  is  no  particle  chemical  reaction, 
nucleation,  condensation,  coagulation  or 
break-up . 

Next,  the  Howarth-Dorodnitsyn  transformation  for 
compressible  boundary  layer  is  introduced: 


a  4 


(12) 


u  =  u  7T-  and  v  =  —  vp  u  a  f 

9  Z  P  9  9 


where  s  and  v  are  the  transformed  x  and  y  coordinates  and  f  is 
a  modified  strearaf unction.  Using  Eqs.  (12,  13)  and  introducing 
nondimensional  temperature  9  (=  T/T^)  ,  Eqs.  (1-3)  become: 


■-H 


f'  -  49 


"+  f9'  =  0 


(16) 
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where  superscript  '  denotes  d/drj. 

The  boundary  conditions  [ Eqs .  (8-9)]  give 

77  =  0:  9  =  9  f  =  f'=Q 

w 

T)->co:  9  =  1  f'=2 


(17) 


For  the  particles,  Eq.  (4),  with  Eqs.  (6-7),  is  transformed 
into : 


|^Y"  +  V(7j)Y*  +  F  ( T| )  Y  =  0  (18) 

where 

V(r 7)  =  f  +  K  9'/9  (19) 

2 

F  (  t?  )  =  K  (99  -  9'  )/92  (20) 

Equations  (15-17)  can  be  solved  first.  The  velocity 
and  temperature  fields  obtained  can  then  be  used  in  Eq.  (18) 
for  the  solution  of  particle  mass  fraction.  The  decoupling  of 
particle  distribution  from  the  flow  field  computation  is  the 
result  of  assumption  of  dilute  particle  concentration.  We  will 
now  concentrate  on  developing  the  analytic  solution  for  Eq. 
(18)  subject  to  boundary  conditions  (10-11).  For  particle  size 
ranging  from  10  nm  to  a  micron,  the  Schmidt  number  is  typically 
high  and  will  range  from  102  to  106  [13].  Therefore,  1/Sc  can 
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be  regarded  as  a  small  parameter  in  Eq.  (18).  The  quantity  V 
in  Eq.  (18)  represents  the  particle  velocity.  As  shown  in  Eq. 
(19),  V  consists  of  a  gas  flow  component  f  and  a  thermophoretic 
component  K Q'/Q.  In  this  problem,  f  starts  with  zero  value  at 
t)  =  0 ,  increases  monotonically  as  v  increases  and  0'  is 
negative  at  t)  =  0,  increases  with  v  and  reaches  zero  at  large 
7).  Thus  at  some  location  above  the  heated  plate,  denoted  by 
77*,  V  is  equal  to  zero.  This  is  the  particle  stagnation  point. 
As  shown  schematically  in  Fig.  1,  on  one  side  of  77*,  V  is 
positive  and  on  the  other  side,  it  is  negative,  i.e.,  V  changes 
sign  at  77  .  Using  1/Sc  as  the  small  expansion  parameter,  Eq. 
(18)  constitutes  a  singularly  perturbed  differential  equation 
with  a  turning  point  [16-19]. 

III.  Analysis 

The  analytic  solution  for  Eq.  (18)  with  boundary 
conditions  (10-11)  will  be  presented  in  this  section  using  the 
matched-asymptotic  expansion  technique  [20-21] . Equation  (18) 
gives  rise  to  two  regions  using  the  perturbation  method:  an 

outer  region  where  only  convection  (including  thermophoresis) 
is  important,  and  an  inner  region  near  17*  where  diffusion  is 
also  important. 

Outer  region: 

Diffusion  is  negligible  and  Eq.  (18)  becomes 


ell 


V(-I?)Y'  +  F(7))Y  =  0 


(21) 


Applying  the  boundary  conditions,  the  solutions  are: 

0  <  n  <  7]#  ,  Y  =  0  (22) 

•  Y  f°°  Ff'n) 

v  >  v  ,  y-  =  expl  (23) 

oo  T) 

In  v  >  V  ,  V  (  tj  )  >  0,  Eq.  (23)  shows  that  if  F  (r/*)  <  0,  Y  (v *)  -> 

0  which  matches  the  value  from  Eq.  (22)  although  the  gradients 

•  * 

dY/d7]  are  not  continuous  at  75  .  If,  however,  F(tj  )  >  0,  Eq. 

(23)  indicates  Y(t)*)/Y  becomes  unbounded.  Appendix  A  sh  iws 
that  the  sign  of  F(t)*)  depends  on  the  value  of  the  product  of 
Prandtl  number  and  thermophoretic  coefficient:  F(rj*)  =  0  if 

KPr  =  1.  When  KPr  >  1,  the  other  solution  is  singular  at  the 

particle  stagnation  point.  This  singularity  can  be  removed  in 
the  inner  solution  where  particle  diffusion  is  included. 

Inner  region 

In  order  to  simplify  the  analytical  treatment  we 

first  consider  a  region  of  thickness  2Av  around  v*  in  which 
V ( tj)  and  F ( 7j )  can  be  linearized.  Eq.  (18)  can  be  rewritten  as: 

Y"  +  (7)  -  7?*)V'  (7J*)Y'  +  |V(7j")  +  (7)  -  7]*)F'  (77*)J  Y  =  0  (24) 
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upon  changing  the  variables  to 


this  becomes 

d2Y  dY 

+  +  ^A  +  b^y  =  0  (25) 
df;  ^ 

with  boundary  conditions 

Y(-l)  =  0  ;  Y  (  +  1)  =  Y^ 

where 

c  =  1/ j^ScA7)2V' (7)*) 

A  =  F  (t?')/V'(7)*) 

B  =  AtjF'(V)/V'(7 )‘)  (26) 

y.,  - 1  exp  [{ .  do,»] 

T7  ♦  A77 

and  £  =  0  corresponds  to  t?  -  tj*.  Near  tj*,  V(tj)  »  0  and  the 
first  term  in  Eq.  (25)  can  only  be  balanced  by  the  last  term  so 
that  an  inner  layer  of  thickness  S-O(v'c)  will  exist.  Outside 
this  inner  layer  and  in  -1  <  (  <  +  1  a  linearized  outer 

solution  will  be  given  by 

Y  =  0  for  -1  *  $  <  0 
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Y  =  Y  ^'Aexp 


[-b(c-i) 


for  0  <  ^  <  +  i 


(27) 


which  can  be  rewritten  in  terms  of  i): 


-A 


Y  =  G  |  T7-77  |  exp  f-  -F 

L  V'  (7)  ) 


t  ,\ 

v-v 


(28) 


where 


G  =  Y  At?V 

♦  i 


and  it  can  be  shown  that  G  is  independent  of  the  choice  of  At? 
within  the  neighborhood  of  tj*  [21].  Note  that  we  must  have 
At?  >  5 ,  so  that  from  the  definition  of  c  in  Eq.  (26)  we  get 


[.  <13 

ScV'(Tj*)J 

Now  we  introduce  a  stretched  variable,  v  =  Z/Vc,  to  transform 
Eq.  (25)  into 


d2  Y 
dv2 


+  v 


dY 

dy 


+ 


0 


(29) 


since  e  is  a  small  parameter  Eq.  (29)  can  be  approximated  by 


d2Y 

dv2 


+ 


+  AY 


0 


(30) 
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The  general  solution  to  Eq .  (30)  is 


y  =  e'u  /4  [c  D  (i>)  +  CD  ( it')  1 

1  A-l  '  2  -A  V  J 


where  the  D's  are  parabolic  cylinder  functions  [23]  and  C  ,  C 

n  12 

constants  of  integration. 

Using  the  matching  principle  described  in  [20,  21] 

the  asymptotic  expansions  of  the  D's  [23]  and  the  boundary 

n 

conditions,  the  following  composite  expansions  are  obtained: 


-1  5  C  -  0  (v  <  V  ) 


Y_(£)  =  Y+iexp(B)  (v'e)"Aexp(-£2/4e) 


m- 


exp  1-A)  Ini  D  £/Vc  +  iAD_A  iC/v^e]  (32a) 

*  -l  J  *  V  J 


0<^5l  (7)  >  7)  )  , 


Y+«-)  =  Y+iexp(B)  (✓c)'Aexp(-?2/4c)  exp [[i-a] rrij t  [$/✓?] 

.  y /  ^  V 

+  iAD_A(iC/v^|  +  Y+iexp(B)  (C)'A[exp(-BO-l]  (32b) 


and  using  the  fact  that 
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D  (iz)  =  r 

/ITT  L 


exp(n7Ti/2)D  (Z)  +  exp 

n- 1 


-run 


(33) 


we  obtain: 


Y  (O  =  Y  exp (B)  (i/e)"A  El^)exp{.?^4c)D  (-$/✓£) 

/ITT 


(34) 


and 


Y+  ( C)  =  Y_(£)  +  Y+iexp(B)  (O 


-A 


exp(-BC)  -  1 


(35) 


this  can  be  rewritten  in  terms  of  tj: 


Y_  (7|)  =  G 

ScV'(7?’)l  i^.(1Ta)  exp 

-SCV'(7 )’)  ) 

L  j  /2TT' 

4 

D 

A-l 


SCV'  (7}  )  (71-7?  ) 


(36) 


Y  (7 1)  =  Y  ( 7) )  +  G(7)-T)VA 


exp 


-F'(7)  ) 

v'(7j*) 


-1 


(37) 


Eqs.  (36-37)  give  the  particle  mass  fraction  profile 
in  the  diffusional  inner  layer  near  the  stagnation  point. 
Together  with  the  outer  solution,  Eqs.  (22-23),  they  give  the 
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complete  distribution  of  particle  mass  fraction  in  the  flow 
field.  To  see  that  the  singularity  at  the  particle  stagnation 
point  is  removed,  we  get  from  Eg.  (36)  or  (37): 


Y(V  )  =  G 


2ScV'  (77 


>] 


A/2 


T(l-A). 
2T ( 1— A/ 2 ) 


From  the  above  expression,  we  see  that  Y  (v  )  will 
[/  Sc  |  .  For  any  finite  value  of  Sc,  Y(v')  will  be 


grow  as 

bounded.  In  the  physical  systems,  Sc  will  be  finite  and  the 
present  analysis  is  valid  for  values  of  Sc  up  to  about  106. 
Beyond  this  value  the  inertia  of  the  particles  may  become 
important . 


IV.  Discussion 

In  this  section,  selected  graphical  presentations 
will  be  given  for  the  analytical  solutions  just  obtained.  This 
is  to  illustrate  the  importance  of  Brownian  diffusion  near  the 
particle  stagnation  point  especially  when  the  product  of 
thermophoretic  coefficient  and  Prandtl  number  is  greater  than 
unity. 


Fig.  2  shows  the  particle  mass  fraction  profiles  near 
the  stagnation  point  (tj  )  as  given  by  Eqs.  (36)  and  (37)  for  Sc 
=  10s,  Pr  =  1  and  KPr  =  .6  and  1.5.  Also  shown  in  the  dashed 
lines  are  the  solutions  with  diffusion  neglected.  For  KPr  <  1, 
we  see  that  particle  diffusion  modifies  the  concentration 
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profile — there  is  a  leakage  through  the  stagnation  point  and 
the  concentration  gradient  drops  smoothly  toward  zero.  This  is 
similar  to  finding  in  Ref.  11.  More  drastic  is  the  case  for 
KPr  >  1  (=  1.5).  The  dashed  curved  in  the  bottom  half  of 
Fig.  2  shows  that  without  considering  diffusion,  the  mass 
fraction  becomes  unbounded  at  the  stagnation  point.  When 
diffusion  is  included,  however,  a  peak  of  Y/Y^  occurs  slightly 
upstream  of  the  stagnation  point  t)  ,  the  particle  concentration 
then  drops  steeply,  passing  through  7)  ,  to  zero.  Thus,  we  see 
that  particle  Brownian  diffusion,  although  having  a  very  small 
diffusion  coefficient,  produces  drastically  different 
concentration  profiles  when  KPr  is  greater  than  unity. 

Fig.  3  shows  the  particle  profile  as  a  function  of  n 
from  the  wall  boundary  to  the  free  stream  for  Pr  =  1,  K  =  1.5 
(Pr  x  K  =  1.5)  and  Sc.  =  104.  In  this  figure,  not  only  the 
inner  solutions  and  also  the  full  range  of  outer  solutions  are 
utilized. 

In  Fig.  4  the  particle  mass  fraction  profiles  near 
their  stagnation  plane  are  shown  for  Sc  ranging  from  102  to  106 
for  Pr  =  1  and  K  =  1.5.  As  Sc  increases  the  peak  sharpens,  the 
peak  value  increases  and  the  peaking  location  moves  closer  to 
n  • 

Fig.  5  shows  the  particle  mass  fraction  profile  for 
different  values  of  K  for  Sc  =  104  and  Pr  =  1.  We  see  that 
with  the  introduction  of  Brownian  diffusion  the  behavior  of  the 
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solution  changes  monotonically  and  smoothly  as  KPr  changes  from 
values  less  than  one  to  more  than  one.  We  also  see  that  as  K 
decreases,  the  particles  stagnation  plane  moves  closer  to  the 
hot  boundary  due  to  the  decrease  in  thermophoret ic  velocity. 

The  analytic  solution  given  by  the  asymptotic 
expansion  has  also  been  compared  with  direct  numerically 
integrated  solutions  of  Eq.  (18).  For  Sc  =  102  to  104  the 
solutions  from  the  two  methods  are  practically 
indistinguishable.  As  expected,  the  asymptotic  solution  breaks 
down  for  moderate  values  of  the  Schmidt  number.  At  Sc  =  10  the 
difference  between  numerical  and  asymptotic  solutions  becomes 
appreciable  and  in  the  analytic  solution  a  slope  discontinuity 
at  Tj  becomes  visible  due  to  the  limit  of  validity  of  the 
expansions  of  the  D^'s.  Direct  numerical  solution  of  Eq.  (18) 
for  Sc  beyond  104  becomes  increasingly  difficult  because  of  the 
demand  of  increasingly  smaller  grid  size. 
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Nomenclature 


A  =  F(t )'}  )  /V'  (v') 

a  stagnation-point  potential  flow  velocity  gradient 

B  -  At?F'  (g*) /V(T)*) 

F  defined  in  Eq .  (20) 

G  defined  in  Eq.  (28) 

f  nondimensionalized  stream  function 

K  thermophoret ic  coefficient 

Pr  Prandtl  number 

Sc  Schmidt  number 

T  absolute  temperature 

V  nondimensional  velocity  Eq .  (19) 

V  particle  mass  fraction,  (mass  of  particle/unit 
volume) / (total  mass  of  gas  and  particle/unit  volume) 

y  distance  away  from  the  heated  wall 

Greek.  Symbols 

V  coefficient  of  dynamic  viscosity 

p  density 

7j  similarity  variable  Eq.  (12) 

9  T/T  ,  nondimensional  Temperature 

00 

5  concentration  inner  layer  thickness 

£  (7 ]  -  tj*)At? 

v  stretched  variable 

c  small  parameter,  see  Eq.  (26) 
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Nomenclature  (continued) 


Subscripts 

t  thermophoret ic 

e  boundary  layer  edge 

oo  free  stream  conditions 

V  <  V 
+  V  >  v' 

S uperscr ipts 

•  particle  stagnation  plane  location 
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Appendix  A:  Sign  of  F(ti  ) 


From  the  definition  of  F(7])  we  have 


F(t)  )  =  K 


'00"-  ©' 


9  J7)  =  7) 


(A.l) 


so  F(t)  )  has  the  same  sign  as  (©©"  -  9'  )  •.  Furthermore, 


71  =  7} 


the  energy  equation  for  9  is 


^0  +  f©  '  =  0 


(A. 2) 


and  at  t)  ,  V(rj  )  =0  therefore  f(v)  =  -  K 


0 


7J  =  7] 


(A. 3) 


Upon  replacing  Eq.  (A. 3)  in  (A. 2)  we  obtain 


00  =  PrK0 '  at  7i=7] 


(A. 4) 


so  that  (90"  -  0'  )  •  =  0'  (7?  )  [PrK  -  1] 


(A. 5) 


and  from  Eqs.  (A. 5)  and  (A.l)  we  see  that 


F (77*)  |  0  if  PrK  |  1 


(A. 6) 
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